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We extend the pure pseudo-power-spectrum formalism proposed recently in the context of the 
Cosmic Microwave Background polarized power spectra estimation by Smith (2006) to incorporate 
cross-spectra computed for multiple maps of the same sky area. We present an implementation 
of such a technique, paying particular attention to a calculation of the relevant window functions 
and mixing (mode-coupling) matrices. We discuss the relevance and treatment of the residual E/B 
leakage for a number of considered sky apodizations as well as compromises and assumptions involved 
in an optimization of the resulting power spectrum uncertainty. In particular, we investigate the 
importance of a pixelization scheme, patch geometry, and sky signal priors used in apodization 
optimization procedures. In addition, we also present results derived for more realistic sky scans 
as motivated by the proposed balloon borne experiment EBEX. We conclude that the presented 
formalism thanks to its speed and efficiency can provide an interesting alternative to the CMB 
polarized power spectra estimators based on the optimal methods at least on angular scales smaller 
than '-^10 degrees. In this regime, we find that it is capable of suppressing the total variance of the 
estimated B-mode spectrum to within a factor of ^ 2 of the variance due to only the sampling and 
noise uncertainty of the i?-modes alone, as derived from the Fisher matrix approach. 
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I. INTRODUCTION 

The reliable characterization and scientific exploitation 
of the polarized Cosmic Microwave Background (CMB) 
Anisotropy signal is one of the main challenges facing the 
CMB research at the present. The challenge is aimed at 
by an entire slew of the experimental efforts which are 
currently on-going [1 , being deployed ^21 131 14l HI 17] or 
planned [8, 9 . The majority, if not all, of these experi- 
ments will effectively observe only part of the sky. This is 
either due to observational/hardware limitations, partic- 
ularly for ground-based and balloon-borne experiments, 
or the presence of foregrounds, or both. The resulting so- 
called cut-sky effects will therefore have to be taken into 
account in the data analysis of most of these experiments. 

The focus of this paper is on the CMB power spectrum 
estimation via the pseudo-spectrum technique p!0| . The 
power spectrum is the most fundamental statistics con- 
sidered in the context of the CMB fluctuations, which 
are thought of as being Gaussian or nearly Gaussian dis- 
tributed. The pseudo spectrum approach provides a com- 
putationally quick and flexible framework for estimating 
the power spectra, which has been extensively and suc- 
cessfully employed in past analyses of the CMB data sets 
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(e.g., [TT1[T2]). However, it has been long recognized [13] 
that a straightforward application of the pseudo spec- 
trum technique to cut-sky polarized CMB maps leads 
to the so-called "£^-to-5" leakage, or power aliasing, as a 
consequence of which the cosmologically important infor- 
mation contained in the CMB 5-modes is overwhelmed 
by the statistical uncertainty of the (much larger) E- 
modes. To date two techniques alleviating the problem 
have been proposed: one correcting the leakage on the 
correlation function level [14] and the other attempting 
to do so directly on the maps [15]. The latter approach 
relies on a suitably chosen sky apodization to remove 
from the map harmonic modes which are neither solely 
E nor B. In a recent paper [16 it has been suggested 
that the apodization can be appropriately optimized to 
bring down the resulting 5-mode spectrum uncertainty 
in-line with the optimal methods, building up a strong 
case in favor of this approach. In this paper we extend 
the original formalism of [15 to incorporate it within the 
framework of the cross-spectra computations and discuss 
a number of implementation issues, related to derivative 
computation and apodization choices in actual applica- 
tions. Later, in this context we discuss the performance 
of the approach. 

We note here that other approaches to the polarized 
power spectrum, which do not suffer from the leakage 
problem, exist but are however hampered either by nu- 
merical efficiency [17] or convergence [18] issues in par- 
ticular for the high resolution observations (but see [19] 



2 



for a potential resolution of this problem) . Consequently, 
the pure pseudo-spectrum approach stands out as an im- 
portant alternative potentially capable of addressing suc- 
cessfully the needs of the future CMB B-mode polariza- 
tion experiments. 

The outline of this paper is as follows. In Sect. [n| 
we start from presenting the standard pseudo-power- 
spectrum technique and the £^/5-leakage problem. We 
also summarize the original approach of Smith (2006) [15] 
and describe its specific variant as proposed in this paper 
and which invokes harmonic domain derivatives and the 
cross-spectrum idea. We conclude that Section with an 
overview of the numerical implementation of the method. 
In Sect, mil we elaborate on the issues related to the 



two approaches are related via the following relations, 



proper sky apodizations. In Sect. [IV| we consider the 
performance of the formalism focusing on the role of the 
apodization and its impact on the quality of the final re- 
sults, as expressed by the level of the remaining residual 
E/B leakage and the uncertainty of the estimated power 
spectra. We summarize the main results of this work in 
Sect. El 

For definiteness, hereafter, we will reserve the term 
mask to a binary, pixelized map of a sky assuming only 
two values 1 and 0, corresponding to an observed or un- 
observed (or just not to be taken into the analysis) sky 
pixel respectively. The pixel weights which can assume 
any nonzero value will be referred to as either apodiza- 
tions or windows. In this language a mask just defines an 
observed sky. We will also refer to aliasing involving dif- 
ferent spin harmonics as leakage, reserving the term alias- 
ing solely to mixing between different modes of the same 
type of the same basis functions. Though not strictly 
precise such differentiations will be useful in the follow- 
ing discussion. 



II. NOTATION AND FORMALISM 

A. Harmonic decomposition 

The CMB (linear) polarization field is completely de- 
scribed by two Stokes parameters, Q and U. Those can 
be combined into two spin-2 and spin- (-2) fields defined 
as, 



P±2 = Q ± iU. 



(1) 



For the full sky, the spin fields can be expressed in the 
harmonic space making use of the spin- weighted spherical 
harmonic basis 



±2 <^£m ±2^^771 



(2) 



Alternatively, the polarization field can be decomposed 
into a gradient, and a curl, 5, part, both of which are 
coordinate independent. In the harmonic space, these 
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From a point of view of physics, the E/B decom- 
position of CMB polarization may seem more natu- 
ral as it is directly linked to the primordial cosmo- 
logical perturbations leading to the presently observed 
CMB anisotropics. In particular, only primordial gravity 
waves, and not density perturbations, can create B-type 
polarization [20^ . As a consequence constraining the am- 
plitude of such modes could allow us to study the char- 
acter of the primordial cosmological perturbations. This 
simple picture is spoilt somewhat by a presence of an 
additional B-mode type contribution generated by the 
gravitational lensing of the E-mode polarization [j21J and 
which needs to be accounted for if the primordial com- 
ponent is to be recovered. 

A computation of polarized angular power spectra 
makes use of spin- weighted spherical harmonics. Those 
functions, a natural harmonic basis for spin-5 fields on 
the sphere, are derived from the standard (or spin-0) 
spherical harmonics by applying the spin-raising (3) or 
spin- lowering (8) operators [20] . 
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where hereafter 5 > 0. In the following, we will always 
assume that s is nonnegative. The main properties of the 
spin- weighted spherical harmonics under different trans- 
formation can be found elsewhere, e.g., [201 l22l [23]. 

Any complex, spin-s field, with s^^^ = -s0 is com- 
pletely characterized by its projection on the spin-s 
spherical harmonics, defining its harmonic representa- 
tion. 



±s<P ±s 



t 



(6) 



Just as is the case for spin-2 polarization fields a spin- 
s field can be also fully described by its E/B harmonic 
decomposition. Adopting a vectorial notation. 



^ = 



Re 
Im 



s4>] 
s4>] 



(7) 



we can define the E and B spherical harmonic basis as. 



= D 
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and 



(9) 




Here we have defined two differential operators Df 
which generalize to arbitrary spin the operators used in 
[13]. The harmonic representation of the field $ in the 
E/B subspace then reads, 



J 



(10) 



where X = E oi B and we have introduced an arbitrary 
(scalar) weight, VK, which we will later need to account 
either for the noise present in the data, or the cut-sky 
effects. The dot denotes a standard componentwise dot 
product. In the special case of CMB polarization, the 
vector $ is simply the polarization vector defined as. 



P = 



U 



It can be easily verified that introducing the polarization 
vector, P into Eq. (10) and setting s = 2 and W = 1 over 



the full sphere, leads to the standard definition of the E 
and B multipoles. 



B. Pseudo power spectra 

The pseudo-Q estimator, e.g. [HI |12j |26] , in a given 
bandpower a with a width A£, is then based on the mul- 
tipole decomposition Eq. ( 10 ) and defined as. 



r)X^ a. 
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where X stands for E or B and i and j distinguish two 
sets of measurements (maps) of the sky for which the har- 
monic coefficients have been computed. Throughout this 
paper, we use i{£ + l)/27r as ^-weighting while comput- 
ing the bandpowers. Though, this may not be the most 
natural choice from the B mode power spectrum per- 
spective, the results presented hereafter are expected to 
be largely independent on it, given the relatively narrow 
^-bins used in this work and the fact that the B-mode 
power spectrum is nearly structure free. For definiteness 
we will be assuming hereafter that the two different sets 
of data always refer to the same sky area, though we will 



allow for different (nonzero) weights assigned to each of 
them. The pseudo-spectrum formalism can be used to 
estimate EE and BB as well as EB spectra. Hereafter, 
we will use a calligraphic typeface, C, for the estimated 
and a roman one, C, for the true underlying CMB power 
spectrum. 

Because of the limited sky coverage, and nonuniform, 
pixel-dependent weights, the above estimator is biased 



and its average over CMB realizations, (C 



involves 



a mixing between different £ modes (or bins) and polar- 
ization states. The latter can be described by a so-called 
mixing kernel which we will denote as . The un- 

biased estimator C^'^.^ is thus obtained by inverting the 
following linear system. 



j^diag j^off 




^E, a 
^B, a 



B, a 



Here N^' 



(12) 

is the noise contribution to the estimated 
pseudo-spectra and vanishes, whenever the noise in the 
two data sets is not correlated, as for example, in a case 
of two data sets produced by two different experiments or 
two uncorelated detectors of a single experiment. In fact, 
we will always assume that N^. .'^^ = 0, whenever i 7^ j, 
and refer to the latter cases as pseudo cross-spectra. We 
will call a pseudo-spectrum computed from a single data 
set, i.e., i = j, a pseudo autospectrum. We note that 
in this case the noise term is usually nonvanishing and 
has to be estimated from the data and/or our knowl- 
edge of the experiment to allow the estimator in Eq. ( 11 ) 



to be correctly unbiased. This emphasizes one of the 
biggest advantages of the cross-spectrum based estima- 
tors, which do not need such a correction. However, the 
mixing kernel, will depend on the window func- 

tions applied to each of the data sets and therefore be 
somewhat more cumbersome to calculate in the cross- 
spectrum case. Moreover, as we will discuss that later 
on, some of the most appealing choices for the window 
function will often require the knowledge of the noise in 
the data. Nevertheless, as the latter will be used only 
for the power spectrum error estimation and not for the 
computations of the spectral estimates, lower precision in 
modeling the noise may be sufficient. In the very least, 
in a case of real data the cross-spectra provide a useful 
and handy diagnostic in the analysis of any data set j27] , 
while for a discussion as the one presented here the cross- 
spectra allow to simplify the analysis by avoiding any ex- 
plicit treatment of the noise bias, while retaining all the 
other aspects essentially unchanged. 

The mixing matrix can be computed as. 
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(13) 



r 



where Xi denotes a position of an z-th pixel. The above 
expression is fully general, as Yx/m can be any function 
used to project the Stokes parameters maps into the har- 
monic E and B subspaces. In particular, for the standard 
spectra the windowing of the data with a window W can 
be taken into account by setting. 



■ Xlm 



(14) 



The specific expressions for relevant for the case 

of standard pseudo spectra can be found elsewhere [12 . 

The matrix is a convolution kernel relating the 

pixel-domain correlation functions of the Q and U Stokes 
parameters to the E and B angular power spectrum (the 
explicit link between two-point correlation functions and 
power spectra for polarization fields can be found in the 
Appendix of [17 ). We point out that in some applica- 
tions considered in the following the window, W ^ and 
therefore also the functions Yx,£m will all depend on the 
bin number, a, (see Eq. (12)). 



In general, the kernel, , will mix different ^- 

modes as well as E and B polarization types. The con- 
sequences of this fact and the conditions, which ensure 
the diagonality of the kernel in the polarization states are 
discussed in the following Sections. 



C. The E/B leakage problem 



The operators defined in Sect. |II A[ valid for any spin-5 
complex field, define two orthogonal subspaces as, 

Df ^ • Df = 0. 

However, the E/B subspaces defined by E/BYsim are 
orthogonal (in the sense of Eq. ([lO|) only on the full sky. 
This means that if the polarization field is known on a 
limited part of the celestial sphere, the reconstruction of 
E and B modes by projecting the measured polarization 
on the E 13^21771 basis is "nonpure", i.e., the estimated 
E or B multipoles receive a contribution from both E 
and B modes. This applies in particular to the pseudo- 
power spectra, C^^ and C^^, as defined in Eq. (11), 



each of which will include a contribution from the other 
mode. Though such a mixing can be removed on average, 
the leaked modes will still contribute to the variance of 
the estimated power spectra. Therefore, if the power 
contained in one of the polarization states is much higher 
than in the other, as is the case of, for example, the CMB 
anisotropics, the error bars estimated for the mode with 
lower power will be drastically exaggerated due to sample 
variance of the leaked contribution. 



The same conclusion can be arrived at by starting from 
Eq. (10) which can be viewed as a harmonic decompo- 



sition of the (lyP) field. Multiplication by the window 
function, W ^ is a convolution not only in the harmonic 
space but also in the E/B subspace, turning thus af^ 
into a mixture of the true B and E CMB multipoles. 

This last observation highlights the fact that the stan- 
dard pseudo-Q technique may be in practice unable to 
provide good estimates of the 5-mode angular power 
spectra in particular from data sets coming from small- 
scale CMB experiments. Moreover, even for satellite ex- 
periments for which angular power spectra can be esti- 
mated from a significant fraction of the sky, the discussed 
effects may be significant enough not to be negligible. 

We can look at the E/B leakage problem also from 
the mixing matrix perspective. The lack of orthogonal- 
ity of E and B type of spherical harmonics results in 
the coupling between the E and B pseudo power spec- 
tra. For this reason the inversion in Eq. ( 12 ) needs to be 



performed simultaneously for both E and B spectra to 
allow for unscrambling the modes with the help of non- 
vanishing off-diagonal EB and BE blocks of the mixing 
matrix. The resulting expression for the power spectrum 
estimate, derived in the case of the 5-mode spectrum 
reads now. 



r>B, a. _ 



i M 

olol' olol 



In the ensemble average sense the above expression is 
unbiased as a result of a subtle cancellation of the E mode 
power present in the pseudo- B and E spectra. Such a 
cancellation does not however apply to the variance of 
the estimator and as a result the variance of the spectra 
of one type will include a contribution from the other. 

In the following we will quantify the E/B leakage using 
both these perspectives, i.e., either looking at an excess 
variance of the estimated B-mode power spectrum or at 
the magnitude of nonvanishing, off-diagonal blocks of the 
mixing matrix. 



D. Pure pseudo-C^ estimators for cross-spectrum 

1. Definition 

Pseudo-C^ estimators of the polarization power spectra 
which do not mix E and B modes can be constructed in 
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projecting the polarization fields on the "pure" E and B 
subspaces [13 . Pure E and B multipoles on a partial sky 
are defined as follows (the vector P now represents the 
polarization signal as measured on the sky) [T5]: 



ECt£m 
BCi£m 



dx[D^{W{x)Yem{x))]^ -P{x), (16) 
I dx [Df (W{x)Y,m{x))] ^ • P(x). (17) 



with W a spin-0 window function satisfying the Dirich- 
let and Neumann conditions on the boundary of the ob- 
served sky region. Such conditions on the window func- 
tion are mandatory for the estimated multipoles to be 
free of a E/B leakage due to partial sky. This can be 
easily shown in the pixel-domain. On applying an inte- 
gration by parts twice to Eqs. ( 16 ) & ( 17 ) and taking into 



account the boundary properties of VK, we can rewrite the 
estimated multipoles as [151 US] i 



BCi£n 



J dx W{x)Yem{x) Df ^ • P(x), (18) 
J dx W{x)Yem{x) ■ P(x). (19) 



Given that the operator (Df ) filters out all the B {E) 
modes, the multipole decomposition defined above indeed 
should be free of any E/B mixing. (We point out that 
the above harmonic decomposition corresponds to the 
decomposition on a partial sky of the so-called x fields 
defined in [H].) A simple way to understand this point 
is to notice that applying to the polarization field 
corresponds to a local filtering of the E modes whereas 
projecting on the b^I im basis corresponds to a global 
filtering. Because of this difference, the latter of these 
two techniques suffers from partial sky effects whereas the 
former one does not. In other words, if the window, 
satisfies the Dirichlet and Neumann boundary conditions, 
the multipoles estimated in Eqs. (16) and (17) are pure 
and the mixing kernel, M^Jf , is a block diagonal matrix 
in the E/B subspace, i.e.. 



XX' 



M^35x,x' 



(20) 



This last condition, intuitively obvious given our previ- 
ous discussion, can be formally derived from Eq. (13) 
defining. 



■ Xa 



= ^^WY,, 



(21) 



A more detailed discussion of the mixing matrix and 
related issues is presented in the next Section and Ap- 
pendix [Aj 



2. Numerical implementation 

a. Pure multipoles. In the actual implementation, 
we do not make a direct use of Eqs. (18) & (19) as this 



would require taking numerical derivatives of noisy maps. 
We rather calculate the pure multipoles using Eqs. (16) 
& (17). We refer hereafter to this way of performing cal- 
culations as the spin-weighted approach. Consequently, 
T>2^^ operators need to be applied directly only to prod- 
ucts of the spherical harmonics times the window func- 
tion. This is an approach also adopted in [1^116|. Our 
numerical implementation proceeds in two steps. First, 
we define two spin- weighted windows. 



Wi = and W2 = S^W, 



and construct three new fields, 



P2 
Pi 
Po 



WiQ^iU), 

wliQ^iU), 

W^iQ^iU). 



(22) 



(23) 



Pg is a spin- 5 field, since = W-s has spin {—s) and 
{Q-\-iU) - spin 2. Then, on the second step, we calculate 
the pure multipoles as linear combinations of the E and 
B multipoles of the three new fields. The latter can be 
derived from Eqs. (16) & (17) with the help of Eqs. (|8| 
& ([9|. For instance, in the case of the 5-modes, the pure 
estimated multipoles read. 




^2,£m+2 



with Bs/m denoting the B-type multipoles of the Pg field 



£m5 



as in Eq. ( 10 ). Numerically the entire computation comes 



down to an efficient calculation of the spin spherical har- 
monic transforms. We comment on those in more detail 



later in this Section. Equation (24) offers another insight 
into the pure estimator formalism. The first term on the 
rhs of this equation corresponds to the standard estima- 
tor, while the two following ones act as counter- terms 
removing the leaked E'-mode. Indeed we have found that 
they are anticorrelated with the first term as expected in 
such an interpretation. 

b. Numerical derivatives. The knowledge of the 
three spin window functions is needed for numerically 
computing the pure multipoles. As will be clear in the 
next Section, there are cases where only the spin-0 win- 
dow is known. Therefore, we need an efficient way to de- 
rive the spin-1 and spin-2 windows from the spin-0 one, 
which in turn requires a numerical calculation of window 
derivatives. Hereafter we will perform this computation 
in the harmonic space. 



(25) 



The final spin-1 and spin-2 windows are then obtained 
by transforming back to the pixel-domain. This tech- 
nique allows us to use simple analytic window functions 
on complicated sky coverage, for which only the spin-0 
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part can be easily computed. It also ensures that the 
conditions given in Eq. (22) are fulfilled. By doing this 
computation in harmonic space, we may introduce some 
ringing due to the pixelization (or gridding) of the sphere. 
Indeed, we have found some minor artefacts of this sort 
when comparing numerically computed windows with the 
respective analytic expressions in cases when the latter 
could be calculated, e.g., in cases with apodizations as in 
Eqs. (30) & (31) and for simple observed sky patches. It 



Diagonal block {N^-^^^=512) 



Off-diagonal block {N^^^^=512) 



appears, however, that such ringing does not seem to af- 
fect our ability to control the E/B leakage at scales where 
this control is mandatory. For instance, for a 20 x 20 de- 
gree square patch pixelized with HEALPix pixels with 
Nside = 512, the residual, leaked modes calculated us- 
ing the numerical spin windows reproduce closely those 
obtained analytically all the way up to ^ ~ 400. For 
higher values of the use of numerical computation leads 
to a slight increase of the residual leakage. We expect 
that to be always the case, independently of the apodiza- 
tion length as long as it is sufficiently greater than the 
pixel size. 

In general the discrepancy between the exact and nu- 
merical spin windows depends on patch geometry. For 
example, the discrepancy is reduced for a spherical cap 
as compared to a square patch. Its magnitude will also 
in general depend on the assumed pixelization and in 
particular may be greatly enhanced if the considered re- 
gion is in a particular position and/or orientation with 
respect to the pixelization scheme. However, whenever 
necessary, such a discrepancy can be reduced or at least 
shifted to higher ^-values by going to a higher resolution. 
We demonstrate some of these effects in Sect. [TVl 

c. Mixing kernel. The mixing kernel can be ex- 
pressed as a sum over the ^-modes in the harmonic space 
of products of the harmonic coefficients of the spin win- 
dow functions and 3j Wigner coefficients [11 . This is 
done by using Eq. (21) in Eq. (13) and subsequently ex- 
panding the window function in harmonic space. We 
present relevant formulae in Appendix [A] (see also [15] 
for an alternative derivation of those kernels). They gen- 
eralize the standard expressions for the ordinary pseudo 
spectra [24l [25] . In our numerical computations we im- 
plement these formulae and perform all the required cal- 
culations directly in the harmonic domain. This is in 
contrast to the previous study of [15 , which strives to 
make computations predominantly in the pixel-domain. 

As we mentioned before for any scalar window func- 
tion we can calculate the harmonic coefficients of the spin 
windows by using the relation in Eq. (25). If we then use 



these in the computation of the mixing kernel, Eqs. (A4) 



& (A5), we find that the off-diagonal blocks of the mix- 



ing matrices vanish within the numerical precision. This 
simply reffects the fact that the pure estimators are in- 
deed "pure" and do not mix E and B modes. This is 
however also not realistic. In practice, we expect that 
some residual E/B leakage will be present due to the 
sky discretization - an effect not corrected by the pure 
pseudo spectra formalism and not characterized properly 
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FIG. 1: Dashed (red) curves show a single column of the mix- 
ing kernel at io = 700 computed using our numerical approach 
as described in the text, for the diagonal, left, and off-diagonal 
blocks, right, and two different resolutions: nside = 512, top 
row, and nside — 1024, bottom. The overlapping solid (black) 
lines depict the EE^ left, and BB^ right, pseudo-spectra com- 
puted from Q and U maps containing only ^-mode power at 
= 700 and apodized with analytic window functions satis- 
fying the boundary condition and constrains in Eq. (22). The 
apodization length is set to 7 degrees. A very good agree- 
ment seen in all the cases validates the proposed numerical 
approach. The off-diagonal elements in the left panels are due 
to the cut-sky mode mixing, while the nonzero values in the 
right panels are due to a residual pixel-induced E/B leakage. 



in the above mentioned calculation of the mixing ma- 
trix. To account for the pixelization effects we propose 
to compute the mixing kernel using coefficients of each 
of the spin windows calculated via numerical spherical 
harmonic transforms (SHT) of the pixel-domain repre- 
sentation, Eqs. (A7) & (A8). We do so even in the cases. 



when the window derivatives are computed numerically 
and thus are directly available in the harmonic domain. 
In such cases, our procedure corresponds simply to tak- 
ing first an inverse SHT of the spin windows as given in 
the harmonic domain. This produces their pixel-domain 
representation, which we transform back to harmonic do- 
main via a direct SHT. Because of sampling effects such 
an operation does not reproduce precisely the input co- 
efficients. The new harmonic coefficients of the windows, 
when used now to calculate the mixing matrix, lead the 
higher values of the off-diagonal block elements, which, 
though still small when compared with the diagonal el- 
ements of the diagonal blocks, are not any more numer- 
ically zero. This can be understood from Eq. (24): be- 



cause of gridding effects during the SHT, the multipoles 
involved in the pure estimators are not computed exactly 
and the precise cancellation of the leaked £^-mode is not 
perfect. 

As those are the pixel-domain representation of the 
windows, which are used in the pure pseudo-spectra com- 
putation, we may expect that the observed change by 
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explicitly taking the SHT of each spin windows while 
computing the mixing kernel, goes in the right direction, 
properly reflecting the correlation pattern between the 
modes of the pure spectra. We can test this expectation 
with the help of numerical experiments. We first note 
that we can compute any of the columns of the mixing 
matrix by computing the pure E and B pseudo-spectra 
of the sky signal containing power only in one selected 
harmonic mode £ = Iq. Examples of such computations 
are shown in Fig [l] There we have allowed for nonzero 
power only in a single io = 700 mode of the E polariza- 
tion. The derived pure E and B pseudo spectra are con- 
trasted with the corresponding columns of the EE and 
BE blocks of M^^/, both computed using the procedure 
outlined before. 

For a square patch assumed in the calculation visual- 
ized in Fig. [ij the overall agreement between the pseudo 
spectra and the column of the kernel is indeed very good. 
For the off-diagonal block some differences can be seen 
only at the high-£ end with the mixing kernel underes- 
timating the residual leakage. The discrepancy appears 
mostly benign as it affects the multipoles for which the 
leaked signal is found to be negligible. 

We thus conclude that the simple way of incorporat- 
ing the pixel effects as proposed here works indeed very 
well everywhere, where the leakage is significant and 
thus requires a sufficiently accurate treatment. This also 
demonstrates that the method of computing the mixing 
matrix advocated here adequately describes the actual 
correlations between the ^-modes of the computed polar- 
ized pseudo-spectra. For more relevant examples of such 
a comparison see Sect. |III A| and Fig. [2j 

We note that in general different window functions can 
be used for the estimation of the power in different £- 
bins. Indeed, as we discuss that later on, it is a common 
occurrence in practical applications in which we usually 
strive for the smallest attainable uncertainty on a final 
power spectrum. In such cases, we calculate the kernel 
row- by-row, i.e. bin- by-bin, using an appropriate window 
assigned to each of them. The kernel computed in this 
way is later used to solve the inverse problem in Eq. (12). 

d. Spherical harmonic transforms. At the core of 
the numerical implementation of the formalism presented 
in the paper is a calculation of the spin spherical har- 
monic transforms. For this purpose we utilize here a pub- 
licly available, Scalable Spherical Harmonic Transforms 
(S^HAT) library [30]. This parallel library allows for an 
efficient computation of arbitrary spin transforms for a 
general class of sky griddings and pixelization in which 
the pixels (grid points) are distributed along the con- 
stant declination lines. These include therefore the most 
frequently used in the CMB work, e.g., HEALPix [31 , 
GLESP [32^, ECP [33] as weh as some others. The 
spin spherical harmonics are computed using the recur- 
rence relations as listed in [34, which are appropriately 
rescaled to avoid numerical under- and overflows. The li- 
brary is designed to treat optimally the cases with only a 
partial sky coverage and is therefore particularly suitable 



for our purposes here. It is parallelized using the Message 
Passing Interface (MPI) and distributes over the pro- 
cessors both the pixel (maps) and harmonic (harmonic 
coefficients) domain objects permitting to achieve high 
resolutions with modest memory resources available per 
single processor. 



III. SKY APODIZATIONS 

Performance and in a result applicability of the pure 
(cross) spectra formalism is clearly dependent on our 
ability to compute a suitable apodization for any specific 
CMB polarization experiment. In this Section we intro- 
duce a number of proposal functions, which could be used 
as apodizations in calculations of the pure (cross) spec- 
tra. The discussed functions range from ones derived as 
a result of optimization procedures and therefore often 
numerically involved, to ones proposed somewhat ad hoc 
but quick and easy to calculate. We will discuss their 
relative merit only briefly here, postponing a more thor- 
ough comparison of their performance to the following 
sections. 



A. Variance-optimized apodizations 

In realistic applications we would like to have an 
apodization function which would ensure that the total 
(i.e., noise plus sample variance) uncertainty of the re- 
covered power spectra is (nearly) minimal. In the case of 
the auto pure pseudo spectra a general ansatz for such 
optimized weighting has been proposed recently by Smith 
and Zaldarriaga [16 . It is derived as a condition on the 
apodization window required to minimize the difference 
between the auto pseudo-spectrum of a map and its opti- 
mal counterpart. Those authors have obtained a follow- 
ing relation for the (nearly) optimal weights. 



Nobs 

E 



1, for all j. 



(26) 



Here the summation is made over the Nobs observed pix- 
els. The matrix C is the covariance matrix of the data 
and P is a geometrical matrix which comes from the 
projection on harmonic space (for example the Legen- 
dre polynomials matrix in the case of temperature) . The 
optimized weight function can be then obtained by a nu- 
merical inversion of the above linear system for each bin 
a. It can be shown [16] that such optimized windows 
converge towards inverse square-noise weighting in the 
noise-dominated case. This is a regime usually reached 
at the high-£ end of the estimated power spectrum. We 
refer the reader for more details about the derivation of 
this relation to [T6] . 

For the computation of the cross-spectrum we need to 
compute an apodization for each of the available maps. 
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The approach of [16] has to be therefore extended to ac- 
commodate such cases. This can be done in a number of 
ways. For instance, using Eq. (26) we can compute an 



Diagonal block (N^ijg=512) 



apodization for each map considering them separately or 
together. Alternately, we could apply a reasoning similar 
to that of [16] but in the context of the cross-spectrum 
to derive a relation analogous to Eq. (26) but applica- 



ble to the cross-spectra. Though this last approach may 
seem the most appropriate for the case at hand it un- 
fortunately results in an expression, Eq. (B14), which 
is computationally rather involved. The two first op- 
tions mentioned above are numerically more tractable, 
though they produce windows, which minimize the dif- 
ference between the pure pseudo-spectrum estimator and 
its optimal counter-part on the auto spectrum level. The 
first of these two approaches uses for that purpose the 
auto spectrum of the map for which the apodization is 
sought, while the second one - the single auto spectrum 
of all the maps considered. Clearly, which one of the two 
methods is more appropriate will depend on the specific 
goals of the entire estimation procedure. We note for ex- 
ample that in the case of two maps with the same noise 
properties the difference between the two optimization 
approaches is limited only to a different overall factor 
determining the level of the noise with respect to the sky 
signal (see Appendix [B| for more details). Hereafter, we 
will utilize the first option and solve Eq. (26) as in the 



pure pseudo auto spectrum case for each of the two maps. 
We will support this choice in the following Sections with 
the help of numerical experiments. 

Equation (26) can be specialized to a particular case 



of polarization when the input data used are composed 
of the Stokes Q and U maps. The geometrical matrix 



>(«) 



needs to be then extended into the spin subspaces 



and nonzero spin windows introduced. Their implicit 
definition consequently reads. 



{a) 



(27) 

where an additional summation over spins of the windows 
is added. The geometrical matrix, P^"), describing the 
projection on spin- weighted spherical harmonics, is now 
explicitly dependent on the spin value, 5, and in general 
is not diagonal in the spin subspace. The vector e is a 
5 Nobs vector whose first Nobs entries are equal to one, 
corresponding to the normalization of the spin-0 compo- 
nent, and all the others are set to zero [TB] . 

The approach can be also looked at as an attempt to 
minimize the total, signal plus noise, power aliased to a 
given bin either from the other adjacent bins or the other 
modes (e.g., E/B leakage). This viewpoint is referred to 
as the variational interpretation in [16 . Consequently, 
in general there is always a different window function 
calculated for each bin as a result of this procedure. This 
has two consequences. First of all, one has to deal with 
a set of windows optimized for each bin. Second, this 
means that one has to define each bin parameters before 




Off-diagonal blocl< {N^^^^=512) 
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FIG. 2: As figure [T] but computed assuming numerically opti- 
mized pixel-domain window functions and a more realistic sky 
patch as could be observed by CMB B-mode b allo on-borne 
experiments, such EBEX, (see Fig. 21 and Sect. IV for more 
details) . 



performing the optimization procedure. 

Below we discuss implementation of this kind of opti- 
mization procedure in both pixel and harmonic domains. 



1. Pixel-domain implementation 



The matrix operations in Eq. (27) should be done 
preferably in a pixel-domain. That allows to consider 
a full potential complexity of the map noise properties 
as well as arbitrary masks with complicated boundaries 
and/or holes due to point-sources extraction, rendering 
such an implementation very flexible. 

In the context of the pure pseudo-spectrum estimators, 
one should in principle ensure that the spin window func- 
tions fulfill the relation in Eq. (22) as well as boundary 



conditions. However, hereafter following [16], we relax 
these constraints in the numerical implementation of the 
procedure. This allows for greater flexibility and there- 
fore enables finding window functions better suited for 
our needs. Nevertheless, this simplification has a major 
consequence: the pure pseudo-Q estimator is not per- 
fectly pure and will unavoidably introduce some E/B 
leakage in addition to the one induced by the pixeliza- 
tion itself. This residual leakage may become problem- 
atic in the low-^ region where sampling variance is the 
major contribution to the error bars. We note however 
that in the high signal-to-noise ratio regime the resulting 
windows are expected to nearly fulfill all the constraints, 
even if those are not explicitly imposed. This is because 
the total variance minimization approach used in J6] to 
derive Eq. (26) implies that the solution for which the 



counter-terms in Eq. (24) nearly cancel the leaked sig- 



nal included in the first term is preferred over those for 
which three terms are independent, which unavoidably 
boosts the variance. In the noise-dominated regime, the 
leaked E/B signal is quickly becoming subdominant and 
therefore it will be the noise variance, not the leaked E 
signal, which will be the major target of the optimiza- 
tion. In this case, the perfect apodization of the window, 
i.e., as required by the pure formalism, may not be the 
best solution, which should gradually, with a decreasing 
signal-to-noise ratio, approach the standard inverse noise 



weighting. 

We point out here that the window function satisfying 
Eq. (26) ensures only that the uncertainty is minimized 
on the pseudo-spectrum not on the power spectrum level. 
In particular, if the mixing matrix is not block diagonal, 
i.e., some residual E/B leakage allowed, the final power 
spectrum uncertainty may not be strictly at its minimum. 
Nevertheless, though some E/B leakage may be indeed 
present, it will be correctly accounted for in our calcu- 
lation of the mixing matrix, and therefore consistently 
treated in the presented method. We show in Fig. |2] that 
the computational approach to the calculation of the mix- 
ing matrix proposed here, and already discussed in the 
case of a simple sky patch and analytic window, fares in- 
deed very well also in the case with the nearly optimized 
window functions. 



The fact that Eq. (27) is solved relaxing the deriva- 



tive constraints between the spin-5 and spin-0 compo- 
nents does not reduce it to a set of three independent 
equations - one for each spin-weighted window. This 
is because the geometrical matrix P-^^ is not diagonal 
in the spin subspace. Consequently all three windows 
have to be solved for simultaneously. In order to do so, 
following [16 we have implemented an iterative Precon- 
ditioned Conjugate Gradient (PCG) solver [35]. The nu- 
merical cost of such an operation is high and on order of 



O {riite 



), if no symmetries are present. A number 



of the iterations, n^ter, is usually as high as a few hun- 
dreds and it is strongly dependent on a dynamic range 
spanned by the pixel noise amplitudes. Nonetheless, for 
small-scale experiments with the help of a proper thresh- 
olding the computational load is clearly within the lim- 
its of available presentday (super) computers. We note 
that for the high signal-to-noise cases imposing explicitly 
the constraints on the window functions could reduce the 
number of iterations required for the procedure to con- 
verge. We have not however attempted to include such 
an optimization in the current version of the software in 
a general case (see however the homogenous noise discus- 
sion in the next Section). 

To illustrate the results of our numerical optimization, 
in Fig. |3] we show three windows computed for a spheri- 
cal cap with homogeneous noise and two different ^-bins. 
The noise level is set to 5.75 /iK-arcmin, the £^-mode 
signal is obtained with the WMAP 5-year cosmological 
parameters [36] and the 5-mode signal includes lensing- 
induced and primordial 5-mode {T / S = 0.05). The ra- 
dius of the cap is roughly 11 degrees, leading to a sig- 
nificant level of E-to-B leakage in the represented bins. 
The resulting window functions are well but not perfectly 
apodized as they do not precisely vanish at the edge of 
the cap. However, as expected, the apodization length, 
defined as a width of the boundary layer, see Eq. (30), 
decreases for bins with higher £. This is in agreement 
with the discussion of [15] concerning the choice of the 
apodization length for homogeneous noise cases. 



spin-g Window 




FIG. 3: From left to right: spin-0, spin-1 & spin-2 spin- 
weighted optimized windows for two bins i G (20, 60) and 
i G (60, 100) are shown in the upper and the lower panels 
respectively. In all cases only the real parts of the spin-1 
and spin-2 windows are displayed. The noise level is set to 
5.75 /xK-arcmin, the ^-mode signal is the WMAP 5-year best 
fit [36] and the B-mode signal includes lensing-induced and 
primordial S— mode {T/S — 0.05). The sky patch is a cap of 
a radius of 11 degree leading thus to a significant amount of E- 
io-B leakage. Note that in the middle (left) panel, the color 
stretch is 3 (6) times bigger for the case with i G (60, 100) 
than for i G (20,60). 



2. Harmonic domain implementation 

If the noise in the pixel-domain is uncorrelated and 
homogeneous over the observed sky patch, the overall 
optimization procedure can be simplified by perform- 
ing the calculations in the harmonic space. That can 



be done most straightforwardly starting from Eqs. (18) 
& (19) rather than Eqs. ([16]) & (17) and therefore results 



in windows for which the boundary conditions, Dirich- 
let + Neumann, are explicitly imposed. 

We begin deriving the average B-mode pseudo- 
spectrum (see Appendix O, 



E 



(28) 



where, 



C} 



= E 



2i' + 1 {£ - 2)\{i' + 2)\ 
■2)!(£'-2)! 



4tt 

Bu2 I 







Here denotes the noise power spectrum, 5^ - the effect 
of an azimuthally symmetric beam. The power spectrum 
C^^// represents the total aliased signal and noise power 
leaking from the £' bin into the £ one. 

The optimized apodization is then given by applying 
the variational principle to wq^ i"m" with three exter- 
nal constrains, which have to be imposed: the window 
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has to be normalized within the observed region and to- 
gether with its first derivative, they have to vanish at the 
contour of the observed region. The computation then 
consists of solving a differential equation with the appro- 
priate Dirichlet + Neumann boundary conditions thus 
leading to the final expression for the optimized window 
given by, 



1 



(29) 



Here M.£"rn" and Vi"m" are the harmonic representation 
of the mask and the contour and A and ji are two con- 
stants given in Appendix [Cl which both depend on the 
signal and noise power. 

In practice the optimization proceeds in two steps. 
First, Eq. (29) is used to provide the scalar component of 



the optimized window. This requires taking the spherical 
harmonic transform of the mask and of its contour and 
calculating the weights Cf^,,. Subsequently, the spin-1 
and spin-2 windows are derived from the scalar one in the 
harmonic space using the numerical tools. This proce- 
dure requires taking the SHT of the binary mask and the 
one of the contour of this mask. Because this last SHT 
is prone to a numerical error (as it corresponds to calcu- 
lating an integral of a pixelized line) , the resulting scalar 
window function is not perfectly apodized. As a conse- 
quence, the window derived in this way does not ensure 
that the estimator is perfectly (or sufficiently) pure and 
may therefore reintroduce E/B leakage. We point out 
that unlike the pixel-domain implementation for which 
the level of leakage is by design smaller than the noise 
level, and therefore under control, the extra-leakage in- 
troduced is here not controlled and thus may or may not 
be smaller than the noise. To cure this problem, we add 
an extra apodization to the scalar window. Wo, before 
computing its first and second derivatives. This extra- 
apodization is obtained by multiplying the spin-0 window 
by the C^-window with a small apodization length (see 
the following Section ). This makes our window optimal 
throughout the entire mask (its shape is driven by the op- 
timization procedure) except very close to the boundary, 
where its shape is driven by the extra apodization. 

The results of the numerical implementation for a 
spherical cap with homogeneous noise are shown in Fig.|4] 
These windows show the same behavior as the windows 
obtained from the PCG with respect to bins and sig- 
nal level. In particular, we find that the apodization 
length decreases for higher bins. The shapes of spin-1 
and spin-2 windows computed in the pixel- and harmonic 
domains are however significantly different. Because of 
the extra-apodization and numerical issues, the windows 
obtained in the harmonic space in a way proposed here 
are not any more an exact solution of the optimization 
problem fulfilling the boundary constraints. Nevertheless 
they provide an easy and quick-to-compute apodization, 
which one may hope will result in spectral variances not 
substantially larger than the one produced by the com- 
putationally heavy pixel-domain optimized windows de- 




FIG. 4: Left to right: spin-0, spin-1, & spin-2, spin- weighted 
optimized windows computed in the harmonic space for i G 
(20, 60) without, top row, and with, middle row, the extra- 
apodization as discussed in the text. We show here only the 
real parts of the spin-1 and spin-2 windows. The color stretch 
of the right panel is 2 times bigger from top to middle row. 
Bottom: The same as middle row but optimized for a bin 
i G (60, 100). In the middle and right panels, the color stretch 
is set to be 4 (2) times bigger, respectively, as compared to 
the corresponding panels of the middle row. In all the shown 
cases the B-mode signal includes lensing and primordial con- 
tribution for T/ S — 0.05 and all the other parameters are set 
as in Fig. [3] 



scribed below. We will test these expectations against 
simulations in Sect. HVl 



3. Sky signal prior and the optimized windows 

The computations of the properly optimized window 
functions requires some prior knowledge about the signal 
expected in the data. In the context of the small-scale 
CMB observations at the low and intermediate i the op- 
timization is mostly between the variance of the E po- 
larization signal leaked to B and the assumed (and pre- 
sumed known) noise level. Hence this is the level of the 
£^-mode polarization, which is of primary importance. As 
polarization experiments measure usually both modes of 
polarization, therefore a high quality internal constraint 
on the £^-mode power should be available for an experi- 
ment aiming at the B-mode detection. At the higher- 
i end a reasonable guess about the B-mode power is 
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^=(60,100), Spherical cap 



FIG. 5: Spin-0, spin-1, & spin-2 optimized window functions, 
left to right, computed for the I G (20, 60) bin for three differ- 
ent assumptions about the B-mode signal: no B-mode (first 
row), lensing-induced B-mode only (second row), lensing- 
induced plus primordial B-mode for T/S = 0.05 (third row). 
In all these three cases, the E-mode signal is obtained using 
the cosmological parameters constrained by WMAP 5-year 
data ^36:. In the fourth row the ^-mode signal amplitude is 
set to a hundredth of that value, i.e., on par with the B-mode 
level, and no ^-mode signal is assumed in the fifth row. In 
these last two cases, the B-mode is the same as in the third 
row. The noise and patch properties are the same as in Figs.jS] 
& [4] The shape of the windows depends weakly on the as- 
sumed B-mode, but clearly more strongly on the ^-mode. For 
visualization purposes, in the fourth and fifth row the color 
stretch is 10 times smaller for the spin-1 and spin-2 windows 
than in the first three rows. 



becoming however increasingly more important. Never- 
theless, the dominant B-mode contribution on subdegree 
scales is due to the gravitational lensing of the ^-mode 
polarization and knowledge of the latter should again be 
sufficient. We elaborate more on these observations in 
the next Section. 
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FIG. 6: Radial profile of the harmonic domain optimized 
spin-0 window function for the noise dominated case (red 
curve) and the signal-dominated case (black curve). For the 
signal-dominated case, the primodial B-mode has been set 
to T/S = 0.5 to emphasize the discrepancy between the two 
cases, all the other parameters are set as in Fig. [5] 



Nonetheless as the prior assumptions have an impact 
on the final results they need to be chosen with care. In 
the case of the pixel-domain optimization we illustrate 
the dependence of the recovered windows on the assumed 
sky signal in Fig. [5] considering five different assumptions 
about the signal, while always keeping the noise level 
fixed at a = 5.75 //K-arcmin. In the first three cases we 
assume the E mode power as in the best-fit WMAP 5- 
year model and change only the B signal. The three cases 
correspond to (i) no B-mode, (ii) only lensing-induced 
B-mode and (iii) both lensing-induced and primordial 
B-mode for T/S = 0.05. In the two remaining cases we 
fix the B-mode power as in the case (iii) and change that 
of the £^-mode to be: (iv) a hundredth of the standard 
value and (v) zero. In all the cases, shown in Fig. [5] 
from top to bottom, the windows have been optimized 
for i G (20, 60) where a high amount of leakage has to be 
corrected for. For assumptions (i) and (ii), the B signal 
is lower than the noise for all angular scales whereas it 
exceeds the noise level in the last three cases. The shape 
of the optimized windows slowly varies with the B-mode 
level meaning that the E/B leakage from the dominant 
E-mode is the main source of extra variance and the op- 
timization procedure focuses on resolving this particular 
source. In the last two cases, the £^-mode is first set to be 
comparable to (case iv) and then lower (case v) than the 
B-mode power, what causes clearly discernible changes 
in the shape of the computed apodizations. In the next 
Section we discuss consequences of the observed changes 
in the window functions on the recovered power spectra 
variance. 

In the case of the harmonic window the optimized 
weight functions depend on the noise and B-mode sig- 
nal prior only, as emphasized by Eq. (29). The B-mode 
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power spectrum appears explicitly in Cf^,, as well as in 
fi and A. Unlike in the pixel-domain case, the £^-mode 
signal is not involved here as by design the harmonic win- 
dow perfectly fulfills the Dirichlet + Neumann boundary 
conditions, ensuring that the £^-mode power leaked to 
B vanishes up to pixel effects. The dependence of such 
a window on the noise and 5-mode prior is depicted on 
Fig. [6j The figure displays the radial profile of the spin-0 
window in the noise-dominated (red curve) or the signal- 
dominated case (black curve). In the noise-dominated 
case, the power spectrum in C^^„ behaves as and our 
numerical result coincides with the analytic expression 
given in Eq. (33) of [16] applied to a spherical cap as the 
survey region, small enough for the fiat sky approxima- 
tion used in [16 to be valid. The apodization length of 
the window in the signal-dominated case is smaller due 



to the fact that in this case C? 



(X ■ 



and thus the power 



spectrum in C^^„ behaves roughly like i' . This conclu- 
sion corroborates the analytic investigation of [16 . 
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FIG. 7: Analytic window functions, spin-0, spin-1 and spin- 
2, left-to-right, respectively, calculated using the (upper 
panels) and (lower panels) functions. Only real parts of 
the spin-1 and 2 windows are shown. 



B. Analytic apodizations 

Analytic window functions are in general not well 
suited for providing an optimal performance in the pres- 
ence of the instrumental noise, given their limited ability 
to adapt to the specific properties of the data. How- 
ever, they are quick to calculate, do not require any 
prior knowledge about the sky signals and can fulfill 
the boundary conditions with high precision. They may 
therefore at least be useful on exploratory, initial stages 
of any CMB data analyzes. Moreover, whenever prior 
information about the sky signals is available, the pa- 
rameters of the analytic window functions can be tuned 
via Monte Carlo (MC) simulations providing a compet- 
itive performance to the optimized windows usually at 
much lower costs. The analytic windows are therefore of 
considerable practical interest and in this Section we will 
discuss two specific proposals of the analytic windows, 
elaborating and complementing the discussion presented 
in [15]. 

First of the considered choices is a straightforward 
generalization of the analytic window function proposed 
in [15] extended to be applicable to an arbitrary sky cov- 
erage. Calling the distance between the z-th observed 
pixel and the boundary, i.e., the smallest angular distance 
form the considered pixel to the contour of the mask, we 
can write the relevant expression as. 



cos 



< 6c 
Si > Sc 



(30) 



and is zero outside the mask. Here, 5c is a width of 
a boundary layer of the window smoothly interpolating 
between the window core and the zero on the outside. 
In the following, we refer to it as an apodization length 
and treat it as an adjustable parameter. Hereafter we 



will refer to the window defined in Eq. (30) as the C^- 
window. One can easily check that if 5i = (i.e., the 
pixel is on the contour of the mask), the window func- 
tion is indeed vanishing, as well as that spin-1 window 
which can be derived from it. Therefore, the window 
defined above indeed fulfills the conditions as required 
by the pure pseudo power spectrum formalism. We note 
however that the related spin-2 window is not apodized 
and potentially therefore could be a source of numerical 
problems in particular, whenever SHT are involved. 

We therefore propose a second analytic formula for the 
spin-0 window from which apodized spin-1 and spin-2 
windows can be derived, reducing in principle the ringing 
due to a sharp spin-2 window. With the same notation 
as above and for every pixel in the mask the window is 
defined as. 



27r 
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1, 



So- 



6i < 5c 
Si > Sc 



(31) 



and vanishes outside the mask. We will refer to this 
choice as the C^-window. 

The three spin windows, for each of the two cases, are 
shown in Fig.jTjin a case of a spherical cap. The apodiza- 
tion length is set to 7 degrees. The spin-0 window indeed 
falls more rapidly to zero in the second proposal, thus 
leading to an apodized spin-2 window. Nevertheless the 
two analytic windows are very similar with only small 
differences concentrated around the edge of the observed 
patch. We will discuss an impact of such differences in 
the window shapes on the E/B leakage in the next Sec- 
tion. 

Two remarks are in order here. First of all, the ana- 
lytic windows can account for the presence of holes in the 
mask for example due to masked out point-sources or un- 
observed pixels. The distance Si just has to be computed 
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FIG. 8: Absolute value of the analytic spin-1 windows (upper 
panels) and spin-2 windows (lower panels) derived via the 
analytic (left line) and numerical computation (middle line). 
The plots in the right panels show the respective differences. 
Note that in the right panels the adopted color scale spans a 
range of values 10~^ (up) and 10~^ (bottom) times smaller 
than in the remaining panels. 



between the i-th pixel and the boundary closest to it (ex- 
ternal or internal). Second, for an arbitrary sky coverage 
(including holes), the spin windows need to be clearly 
derived numerically, however, in some simple cases such 
as that of a spherical cap or square patches the spin-1 
and spin-2 windows can be computed analytically from 
the spin-0 window. We can, and will, use that latter ob- 
servation to quantify some numerical effects, involved in 
the discussed formalism. For instance we will investigate 
the effects of the numerical derivatives on the shape of 
the spin- weighted windows, as shown for a spherical cap 
case in Fig. [8) 



IV. TESTS AND APPLICATIONS 

Performance of the technique discussed here depends 
on many factors, ranging from the pixelization effects 
to assumed knowledge of the sky signals, and can be 
looked at from different perspectives. In this Section we 
discuss some of the most relevant effects. As the pre- 
sented estimator first aims at solving the E/B leakage, 
we start this Section with a discussion of a quality of the 
E/B leakage control and follow on with a study of the 
estimator efficiency expressed in terms of the 5-mode 
power spectrum variance levels, which can be achieved 
using this approach in different circumstances. Unless 
it is explicitly mentioned otherwise, we will work with 
the HEALPix pixelization scheme adopting the resolu- 
tion Nside = 512. The input £^-mode signal is that of 
the cosmological model with parameters as constrained 
by the WMAP 5-year data [36 . The 5-mode signal in- 
cludes lensing and primordial 5-mode with T/S = 0.05. 



For the case of homogeneous noise the noise level is set 
to 5.75 /iK-arcmin, which corresponds to a typical level 
expected for the small-scale B-mode experiments. 

To study the dependence on the patch geometry, we 
adopt five different patch shapes: a spherical cap, a 
square and three rectangular patches with a different 
elongation. Their geometrical properties are summa- 
rized in Table [ij All the patches cover the same sky area 
(roughly 1%) but the length of their contour increases 
progressively from the spherical cap (C) to the most elon- 
gated rectangle (R3). 
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(R3) 
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49.8 


2.01 


CAP 


(C) 
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11.3 deg. 


1.23 



TABLE I: Geometrical properties of the five patches used in 
the examples considered here: is the zenithal aperture, cp the 
azimuthal aperture and the perimeter is the total geometrical 
length of the contour of the patch. The five patches cover 
the same sky area roughly equal to 1% of the entire celestial 
sphere. 



A. Mode mixing kernel 

The mode mixing kernel provides the most direct mea- 
sure of the E/B leakage. The root-mean-square (rms) of 
the fraction of the power contained in one of the two po- 
larized modes which is leaked to the other mode is char- 
acterized by the magnitude of the off-diagonal elements 
of the mixing kernel. For the pure estimators in princi- 
ple those should be vanishing if the proper apodization 
window is used. However, in practice, the apodization 
is only one of the factors on which the level of the leak- 
age depends and numerical effects due to approximations 
and simplifications involved in the calculations of the es- 
timated spectra give often rise to the power leakages. 
In this Section we study dependence of the off-diagonal 
blocks of the mixing kernels on the specific choices of the 
apodization windows, pixelization scheme, pixel size, sky 
location and shape of observed sky areas, as well as er- 
rors due to the approximations involved in the numerical 
derivative calculations. 

We show an example of the mixing kernels for both 
standard and pure pseudo-spectra in Fig. [9) In the stan- 
dard case we note a strongly band diagonal character 
of both, diagonal and off-diagonal, blocks of the mixing 
kernel with a magnitude of the diagonal elements of the 
off-diagonal block comprising a significant fraction of the 
corresponding elements in the diagonal block. This is a 
clear indication of the leakage. The off-diagonal block el- 
ements are manifestly suppressed in the case of the pure 
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This leakage appears therefore mostly due to the cut-sky 
effects, which are, by design, not fully removed in the 
case of the optimized windows which do not fulfill the 
proper boundary requirements, and are only suppressed 
to the levels below that of the other uncertainties. 



(r = 80, Cap) 
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FIG. 9: Upper row: Diagonal, BB, (Wt) cind off-diagonal, 
EB, (right) blocks of the mixing kernels for the standard es- 
timations. Lower row: Same as above but for the pure es- 
timator computed using the analytic C^-window. The color 
scale is logarithmic and the same in all the panels. The verti- 
cal axes correspond to the power spectrum multipoles, while 
the horizontal ones to the pseudo-spectrum, Eq. (12). 



estimator, as expected given that in the displayed case we 
have used a properly apodized, analytic C^-window. Also 
as expected the diagonal block remains strongly diago- 
nal dominated. Nevertheless there is a substantial level of 
mixing between the low-^ multipoles of the power spectra 
with essentially all multipoles of the corresponding pure 
pseudo- spectra seen for both the blocks. This extra mix- 
ing arises due to the presence of the counterterms in the 
definition of the pure modes, Eq. (24). As a consequence 



a special treatment on the map level aiming at reducing 
the power contained in the modes comparable and larger 
than the observed sky patch may be required prior to the 
power spectrum estimation. 
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FIG. 10: Off-diagonal block of the mixing kernel M^^/ for two 
different classes of windows and for different HEALPix pixel 
sizes: Nside = 1024 (black-solid curves), Nside = 512 (blue- 
dashed curves) and Nside = 256 (red-dashed-dotted curves). 
The three lower curves are for the -window and the three, 
overlapping upper curves - for the optimized window com- 
puted in the pixel-domain. The apodization length has been 
set equal to 9 degrees for the C^-window and the optimized 
windows has been optimized for i G (60, 100). For the ana- 
lytic windows, the residual leaked ^-mode is induced by pixel 
effects and therefore decreases for the higher resolution. For 
the optimized window, the residual leakage is mainly due to 
a sky cut and as such does not depend strongly on the pixel 



Pixel size 



2. Window choice 



As we already hinted at before (Sect. IIP 2| ), pixeliza- 
tion of the sphere is a source of residual E/B leak- 
age, which will in general depend on the pixel sizes 
and shapes. The dependence of such a leakage on the 
HEALPix pixel size for two types of window functions is 



displayed in Fig. 10 Each of the three lower curves shows 
a single column of the off-diagonal block of the mixing 
kernel, ^ computed using the C^-window and cor- 

responding to, from top to bottom, Ngide = 256, 512 & 
1024. The residual, leaked £^-mode decreases for higher 
resolution confirming the pixel origin of this leakage. The 
three upper curves in Fig. [lO] show the case of the op- 
timized window computed in the pixel-domain. These 
windows have been optimized for the range of I from 60 
to 100 and we kept the noise per pixel area constant. In 
all three cases we find a substantial level of the residual 
leakage, which is independent on the adopted pixel size. 



A comparison of the residual levels of the E-to-B leak- 
age for the different sky apodizations as considered here 
is presented in Fig. [TT] We plot there a single column 



of the off-diagonal block of the mixing kernel for = 80 
and computed for the four types of windows defined in 
Sect. |III| and the spherical cap area as defined in Ta- 
ble [Tj As expected the optimized window computed in 
the pixel-domain leads to the highest amount of leakage. 
The harmonic space window fares better in part thanks to 
the boundary conditions which are explicitly imposed on 
it and in part due to the extra-apodization. For the two 
types of optimized windows, the level of leakage varies 
slightly depending on a prior assumed for the B-mode 
power. 

Both the analytic windows satisfy the Dirichlet and 
Neumann conditions and therefore lead only to a residual 
amount of E-to-B leakage. Nevertheless, the C^-window 
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FIG. 11: A = 80 column of an off-diagonal block of the mixing kernel M^J/ for the four classes of apodized windows. Left 
panel : C^-analytic window for three different values of the apodization length : 6c = 3 deg. (black-solid curve), = 5 deg. 
(red-dashed curve) and Sc = S deg. (blue-dashed-dotted curve). Middle panel : Same as left panel but for the C^-window. 
Right panel : optimized windows computed in the pixel-domain (black curves) and in the harmonic domain (red curves). Solid 
curves are for optimized windows derived without 5-mode and dashed curves for optimized windows with lensing-induced and 
primordial B-mode for T/S = 0.05. 



performance is superior to that of the C^-window, pro- 
ducing the level of leaked £^-mode 10~^ to 10~^ smaller. 
This is because the level of the leakage is reduced by in- 
creasing the apodization length. For the analytic win- 
dows, more apodization is achieved either by increas- 
ing the apodization length or by additionally forcing the 
spin-2 window to vanish at the edges of the survey. For 
the optimized windows, as discussed before the ampli- 
tude of the adopted B-mode prior leads to an increase of 
the effective apodization length, thus affecting the leak- 
age level. Because of pixelization, the windows are dis- 
cretized and the Dirichlet and Neumann boundary con- 
ditions are inevitably broken. By increasing the apodiza- 
tion, we reduce the size of the discrete jump in the win- 
dow values between two neighbouring pixels at the patch 
boundary, thus mitigating the extent to which that af- 
fects the boundary conditions. 
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FIG. 12: A column of a diagonal (upper curves) and off- 
diagonal (lower curves) block of the mixing kernel using an- 
alytically (black-solid curves) and numerically (red-dashed 
curves) computed spin-1 and spin-2 windows. The dashed 
and solid curves overlap over most of the angular scale. Left 
panel shows a column calculated for i' = 80 and right - for 
f — 500. The numerical computation of the spin-1 and spin-2 
windows leads to a slight increase of the leakage in the high 
i part seen as a difference in the off-diagonal block curves. 



3. Derivatives 

For arbitrary patch geometries and in the case of the 
analytic apodizations the spin-1 and spin-2 windows may 
have to be computed from its spin-0 counterpart numer- 
ically and therefore being often only approximate. In 
this paper, as described before, we perform the numerical 
derivation in the harmonic domain. The effect of using 
the approximate rather than exact nonzero spin windows 



is depicted in Fig. 12, where again a single column of the 
mixing kernel is shown for the spherical cap computed 
using either the analytic formula (in black) or numerical 
computations (in red) for spin-1 and spin-2 component. 
This shows that calculating the derivatives numerically 
does not affect the diagonal block elements and leads only 
to a slight increase of the E-to-B leakage, as expressed 
by a small change in the off- diagonal block element mag- 
nitudes. This latter effect is indeed very small and only 
appears at the smallest scales, confirming thus that our 
prescription for the derivative computation permits us to 



retain very good control of the leakage, especially in the 
low-^ part, where it is of primary importance. 



4- Patch geometry 



The remaining leakage as a function of the geometry 



of the observed sky area is depicted in Fig. 13 for the 
two analytic windows and for two values of the apodiza- 
tion length. For an apodization length of 3 degree the 
level of leaked £^-mode is the same in all the considered 
cases, demonstrating that we can achieve a good control 
of the leakage for any elongation. For (5c = 5 degree the 
leakage for the most elongated patch, i.e., the R3 geome- 
try, is greatly enhanced. This is because for such a patch, 
apodization lengths greater than 4 degree exceed the half 
of its smallest size causing the spin-0 window, as defined 
in Eqs. (30) or (31 ), to be strictly speaking nondifferen- 



tiable. Nevertheless the spin-1 and spin-2 windows can 
still be formally introduced, using the fact that one-sided 
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derivatives exist. That however unavoidably results in 
the spin-1 and spin-2 windows, which are not always con- 
tinuous, leading the leakage excess. This emphasizes the 
importance of ensuring that all three spin-5 windows are 
continuous throughout the entire observed region - a con- 
dition, which on occasion may require either trimming 
some irregular edges of realistic sky patches as observed 
by actual experiments or applying an extra smoothing to 
the spin-0 window. Nevertheless, whenever this condi- 
tion is fulfilled the leakage is largely independent on the 
observed patch shape, and with the help of the presented 
formalism can be, at least in principle, removed down to 
the level as defined by the pixelization effects. 
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FIG. 13: A column of the off-diagonal block of the mixing 
kernel M^J/ for t — 80 and for five different patch geometry: 
spherical cap (black-solid curve), square (blue-solid curve), 
rectangle Rl (red-dashed curve), rectangle R2 (black-dashed- 
dotted curve) and rectangle R3 (blue-dashed curve). The 
left panels show the results for the -window and the right 
panel for the C^-window. Upper panels are for 5c — ^ deg. 
and lower panels for = 3 deg. For the (R3) geometry, an 
apodization length of 5 deg. exceeds the half of the smallest 
size of the patch. As a consequence, the spin-1 and spin-2 
windows are not continuous right in the middle of the patch, 
leading to a significant increase of the pixel-induced leaked 
E-uiode (see blue-dashed curves of the upper panels). For all 
the other cases, the curves overlap proving the good control 
of the E/B leakge irrespectively of the geometry. 

The case of optimized windows computed in the pixel- 



domain is displayed on Fig. 14 where the level of leakage 
is shown for the windows optimized in the ^ G (60, 100) 
bin. The noise and signal prior is the same for all the 
patches. As for the analytical sky apodizations, the level 
of leakage only marginally depends on the patch geome- 
try. This means that the optimization procedure manage 
to lower the level of leaked E mode below the noise irre- 
spectively of the shape of the boundaries (we remind that 
the five surveys cover the same sky area, which roughly 
leads to the same amount of noise in £-space). We note 
that as expected the level of leakage allowed by the opti- 
mized windows is a few orders of magnitude higher than 
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FIG. 14: A column of the off-diagonal block of the mixing 
kernel M^J/ for i' = 80 and calculated for windows opti- 
mized in a pixel-domain for a range i G (60, 100) and the 
five different patch geometries (the five curves overlap): 
spherical cap (black-solid curve), square (blue-solid curve), 
rectangle Rl (red-dashed curve), rectangle R2 (black-dashed- 
dotted curve) and rectangle R3 (blue-dashed curve). The axis 
ranges coincide with those in Fig. [13] 



for the analytic windows. 



5. Dependence on a sky position/pixelization 

Numerical computations of derivatives of spin-0 win- 
dows will depend in general on the sky location and on 
the adopted pixelization scheme, even if the windows are 
defined analytically. For instance, in the case of a spher- 
ical cap centered at the north pole, the precision of a 
calculation of the azimuthal part of the SHT quickly de- 
teriorates in the HEALPix scheme because of the de- 
creasing number of constant declination pixels present in 
the polar areas. As a consequence, the spin-1 and spin-2 
windows computed in the harmonic space can, sometimes 
drastically, depart from their true behavior often exag- 
gerating the level of the leakage. Such a problem can be 
cured either by rotating the cap close to the equatorial 
plane, as often plausible in cases with a small sky cover- 
age, or by using a different pixelization/gridding scheme 
such as GLCP (i.e., Gauss-Legendre Celestial Pixeliza- 
tion - ECP-like but with the constant declination rows 
distributed as the zeros of the Legendre polynomial of 
some order as supported by the S^HAT software [30]), 
for which the number of azimuthal pixels (grid-points) 
remains the same for each constant declination row, and 
which provide a nearly exact quadrature over the polar 
angle. We provide examples of the potential effects in 
We note that the choice of the pixelization used 



Fig. 15 



for the derivative computations can be different from that 
used for the analyzed maps, though care needs to be 
taken while performing interpolations required in such 
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FIG. 15: An impact of the sky location and the pixeUzation scheme on the diagonal block (upper panels) and off-diagonal 
block (lower panels) of the mixing kernel M^^/. The different curves partially overlap for the diagonal block. Left panels : 
Computation in the HEALPix scheme for a cap centered at the equatorial plane using analytic spin-1 and spin-2 windows 
(black-solid curve) or numerical spin-1 and spin-2 windows (red-dashed curve) and for a cap centered at the north pole using 
analytic spin-1 and spin-2 windows (blue-dashed-dotted curve) or numerical spin-1 and spin-2 windows (green-dashed-dotted 
curve). Middle panels : Same as the left panels for ECP pixelization scheme. Right panels : Same as the left panels for GLCP 
pixelization scheme. The window function used here is the C^-window with the apodization length set 6c = 8 degrees. In all 
three cases the pixelization parameters are chosen in such a way that there is the identical number of the constant declination 
rings within the the sky cap for each of them. For the HEALPix scheme, the bad sampling of the sphere close to the north 
pole (particularly along the azimuthal direction) leads to strong pixelization effects becoming dominant even for the diagonal 
block. A similar problem is also seen in the ECP case, but is resolved in the case of the GLCP scheme, which properly samples 
the sphere all the way to the poles in both polar and azimuthal directions. The effect can be also alleviated by shifting the cap 
closer to the equatorial plane. 



a case. In addition, similar effects may also play a role 
whenever any SHT is to be done, as for example, in the 
computation of the pure multipoles, and therefore choice 
of the pixelization and/or its orientation are of primary 
importance for techniques as the one discussed here. 



B. Power spectra uncertainty 

The errors of the estimated B-mode power spectra in- 
clude contributions due to both the signal and the noise. 
For any ^-bin the signal sample variance is due to the 
variance of the B-mode signal modes from that bin but 
also from other multipoles aliased to it. It is also due 
to power of the leaked £^-mode into 5, i.e., the E/B 
leakage. The window functions discussed before lead to 
different trade-offs between these different sources of the 
uncertainties. In this Section, we will assess the power 
spectra uncertainty using MC simulations comparing the 
performances of the pure pseudo cross power spectrum 
approach for different choices of the window functions. 
We will also consider the effects due to patch geometry 



and signal priors, complementing the discussion in the 
previous Section. 

Hereafter, we estimate the variance as the standard 
deviation of 1,000 MC simulations. Unless specified ex- 
plicitly otherwise we will employ the same set of param- 
eters as before, i.e., we will work with the HEALPix 
scheme with the resolution Ngide = 512 and assuming 
the WMAP 5-year data cosmological model [36 for the 
input £^-mode signal. The B-mode signal includes lensing 
and primordial B-mode with T/S = 0.05. We first focus 
on the homogeneous noise case with a level of 5.75 /iK- 
arcmin. The estimated power spectrum are binned with 
the bins of a width of Ai = 40 and with the lowest-^ 
bin starting at ^ = 20. Except for the cases investigating 
effects of the patch geometry, the simulated maps cover 
a spherical cap (C) area centered at the equator as de- 
scribed in Table [H 

As a benchmark for the derived numerical results we 
will use the simplified, theoretical variance derived for a 
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FIG. 16: Average B-mode spectra (solid curves) and their 
sample variances (dashed curves) obtained using standard, in 
blue, or pure estimators, in red, and computed from 1000 
signal-only MC simulations done on a spherical cap cover- 
ing 1% of the sky. Also shown are the input power spec- 
trum (black-solid) and the Fisher estimate of the variance 
(black-dashed). The extra- variance for the standard estima- 
tors comes from the ^-mode power leaked into B, while the 
sample variance for the pure estimators is close to optimal. 



cross-spectra from the Fisher consideration [TT] [29] , 
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2/ 
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B 2 
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Bf 



(32) 



where a is the noise power spectrum and fsky will be 
hereafter taken to be an actual covered sky area. As such 
this uncertainty provides a lower limit on any power spec- 
trum estimation error, including the one derived using 
the pure estimator discussed here, and therefore a useful 
reference to gauge the performance of the pure method 
against. We emphasize that we do not expect that our 
estimator will be ever able to reach this limit even in the 
cases when no E / B leakage is allowed for, i.e., when Cf 
is set to 0. In fact a tighter lower limit could be read- 
ily obtained, if an effective, and thus usually smaller, 
rather than actual sky area is used as fsky in Eq. (32) 



accounting for the presence of the apodizations in the 
pure formalism. However, even then the excess variance 
is usually seen owing to the presence of the counter terms 
in Eq. (I24|). 



We illustrate some of these effects in Fig. 16, where 
we present the results of the signal-only MC simulations 
for both standard and pure pseudo-spectra approaches 
and contrast them with the semianalytic predictions as 
calculated using the Fisher approach, Eq. (32). 



against the variance derived using the optimized windows 
computed in the pixel-domain (red curve) and its Fisher 
estimate (black-dashed curve). As a reference, the input 
theoretical 5-mode used in the simulations is plotted as 
black-solid lines. Given the adopted here sky coverage 
(and hence bin-width) as well as the noise level, both in 
a ball-park of the anticipated experiments, we find that 
in the range of ^ from 20 to 1000 the variance is always 
few times smaller than the signal itself thus in principle 
allowing for a good reconstruction of the 5-mode power 
spectrum. 

The apodization length has clearly an effect on the 
performance of the method. For the C^-windows, the 
lowest variance is reached using a 5 degree apodization 
in the first bin and 2 degrees in all the others. For the 
C^-windows, 4 and 3 degrees should be used in the first 
and second bins respectively, whereas the 2 degree win- 
dow gives the lowest variance for all the other bins. This 
is due to the fact that noise in the subsequent bins is big- 
ger for the higher i values and thus though the leakage 
level somewhat increases for the smaller apodizations, 
boosting up the variance. This is compensated by the 
simultaneous decrease of the noise variance. At least for 
the homogeneous noise the best optimization length can 
be clearly tuned for every analytic window and each bin 
if some prior information about the sky power is avail- 
able. This can be done determining via MC simulations 
the total estimator variance for a range of considered 
apodizations and selecting the one leading to the mini- 
mal scatter. 

The variance derived using the optimized windows 
computed either in the pixel-domain (red curve) or in 
harmonic space (blue curve) is displayed in Fig. 18 The 
two windows lead to a similar level of variance, which is 
also on par with the Fisher forecast at least for I > 200. 
Clearly, in the cases studied here the use of pure pseudo- 
spectrum technique yields results close to the best achiev- 
able for angular scales smaller than a degree. Moreover, 
the derived variance is indeed at least no larger than 
the one obtained with the analytical windows as can be 
seen in Fig. fTTl The very good agreement between the 



harmonic space computation and the pixel-domain com- 
putation of the optimized windows shows that the extra 
apodization used during the harmonic computation does 
not cause any significant departure from optimality. This 
justifies the use of the harmonic window in the homoge- 
neous, or nearly so, noise cases as it is much faster, and 
better controls the E/B leakage across the entire range 
of i and does it without any apparent loss of performance 
for the final precision of the estimated power spectrum. 



1. Effect of the window function 

The variance of the estimated B-mode power spectra 
for a number of the apodization lengths and the two an- 



alytic proposals is displayed in Fig. 17 and contrasted 



2. Effect of the signal prior 

As we have emphasized before the optimization pro- 
cedure requires knowledge of both the noise and signal 
properties, including the B-mode. In most of the cases 
of the forthcoming CMB polarization experiments noise 
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FIG. 17: Variance of the estimated B-mode power spectrum using the C"^-window, left panel, and the C^-window, right panel, 
for different apodization lengths: from 2 (blue) to 9 (yellow). The red curve shows the variance obtained using the optimized 
windows computed in the pixel-domain and optimized bin-by-bin. The observed survey is the spherical cap (C) and the 
cosmological model and the level of homogeneous noise are described in the text. The input theoretical B-mode is shown with 
a black-solid line and a Fisher estimate of the variance with a black-dashed curve. 
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FIG. 18: Variance of the estimated B-mode power spectrum 
using the optimized windows computed in the pixel (blue) 
and harmonic space (red curve). The black-dashed curve 
shows the Fisher estimate of the variance and the black-solid 
line - the input theoretical B-mode. The same performance 
is achieved using the two implementations of the optimized 
window and is close to optimal for i > 200. 



and £^-mode properties will be probably sufficiently con- 
strained for high and intermediate i modes, so the major 
source of uncertainties will result from our ignorance of 
the primordial B-mode and \ow-£ £^-mode power. 

To assess the influence of the B-mode prior on the 
performance of the optimized windows, we have per- 
formed simulations using lensing-induced and primordial 
B-mode {T/S = 0.05) as input signal but estimating the 
power spectrum using three different sets of windows op- 
timized for different B-mode signal : (1) no B-mode, 
(2) lensing-induced B-mode and (3) lensing-induced plus 
primordial B-mode for T/S = 0.05. As the two first as- 
sumptions differ from the B-mode used to synthesize the 
map in the simulations, we expect that the optimized 



windows derived under those hypotheses could occur to 
be suboptimal. The results of such simulations are shown 
in Fig. 19, left panel, where the signal variance (solid 



curves) and the noise variance (dashed curves) for the 
three sets of optimized windows are displayed. The black 
curves correspond to their Fisher estimates. It clearly 
shows that the prior on the B-mode used in the optimiza- 
tion procedure does not significantly affect the perfor- 
mance of the pure pseudo-spectrum estimation, as long 
as the postulated power of the B-mode remains much 
smaller than that of assumed in the E-mode prior. This 
last condition will be always fulfilled taking into account 
the current upper limits on the B-mode at least as long 
as CMB is the only sky signal considered. As a conse- 
quence, a conservative, but sufficient, i.e., nearly opti- 
mal approach would consist in assuming lensing-induced 
B-mode for optimization, as its level and shape will be 
largely known from the modeling of the £^-mode spec- 
trum. 

Though in general we expect to have a sufficiently 
well constrained £^-mode power spectrum for the future 
polarization-sensitive experiments, that may not be al- 
ways so at least as far as the very lowest multipoles are 
considered. To test the impact of those on the vari- 
ance of the recovered power spectra we have computed 
the optimized windows with a fixed B-mode power spec- 
trum, including lensing-induced and primordial B-mode 
for T/S = 0.05, but using three different E-mode pri- 
ors. The three E-mode power spectra used in this test 
were calculated for three different values of the reioniza- 
tion optical depth r: (1) r = 0, (2) r = 0.2 and (3) 
r = TwMAP = 0.087 as constrained by the WMAP 5- 
year results. These three power spectra have been used 
to first compute the corresponding optimized windows 
and later to perform MC simulations in order to obtain 
the variance estimates. The 'true' sky used to produce 
the maps to-be-analyzed here assumed r = twmap and 
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FIG. 19: Left panel: Noise (dashed curves) and sampling (solid) variance of the estimated B-mode power spectrum using 
pixel-domain optimized windows calculated for three different assumptions about the B-mode signal: no B-mode (blue curves), 
only lensing-induced B-mode (red curves) and lensing-induced plus primordial B-mode for T/S = 0.05 (yellow curves). The 
Fisher estimates of the noise and signal variances are shown in black. This is the 3rd case model, which has been used to 
generate the signal maps analyzed in all three cases. Though the assumed sky priors for B-modes in two first cases are wrong, 
that does not seem to have a strong impact on either the signal or the noise variance. Right panel: Same as left panel but using 
windows optimized for three different assumption of the B-mode signal via a varying reionization optical depth: r = (red 
curves), r = 0.2 (yellow curves) and r = twmap = 0.087 (blue curves). This last model was also used to generate the signal 
maps analyzed in all three cases. As before no strong impact of the priors on either the recovered signal or the noise variance 
is seen in any of the considered cases. 



thus coincided with the third prior as considered above. 
The resulting noise (dashed-colored curves) and sample 
(solid-colored curves) for the three prior are displayed on 
Fig. 19, right panel, alongside the Fisher estimate (black- 



dashed curve) and the input Cf (black-solid curve). We 
find that the variance is only marginally affected by the 
B'-prior, even at large angular scales where a varying op- 
tical depth significantly changes the shape and amplitude 
of the B-mode. We note that our knowledge of the E- 
mode spectrum will be soon much more precise than the 
uncertainty permitted in the test described here, provid- 
ing therefore an external prior certainly sufficiently ac- 
curate for the purpose of the method discussed in this 
paper. 
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3. Effect of the patch geometry 

From the previous discussion, we have already learned 
that the shape of the observed sky area does not affect 
our ability to control the E/B leakage with the help of 
a proper apodization and pure formalism. We can there- 
fore anticipate that the dependence of the variance of 
the estimated B-mode power spectrum on the patch ge- 
ometry will follow the same pattern as for the standard 
pseudo-spectrum estimator. 

We show the numerical results in Fig.|20j alongside the 
Fisher estimate, which is the same for all considered sky 
patches as they all have the same sky area (Table The 
power spectra have been estimated using the -window 
function with an apodization length optimized for each 
bin and every patch separately. Except for the REC3 



FIG. 20: Variance of the estimated B-mode power spectrum 
with respect to the patch geometry using the C^-window with 
an apodization length optimized via MC simulations for each 
bin and each considered patch. The black-dashed curve shows 
the Fisher estimate of the variance and the black-solid line - 
the theoretical input B-mode spectrum. For I > 200, the 
estimation is close to optimal and a good reconstruction of 
the B-mode can be obtained. For smaller multipoles, the 
performance deteriorates for more elongated patches though 
a variance smaller than the expected power level is recovered 
in most of the studied cases. 



patch, for which the variance of the estimated power 
spectrum exceeds the power level in the first bin, the 
variance is small enough to ensure a good reconstruction 
of the B-mode power spectrum. Moreover, for all the 
patches, it is roughly the same and on the level com- 
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inversion of the involved matrix is clearly too costly to 
be performed. For this reason in the following examples 
together with the full observed sky patch we consider 
also its well-observed part as shown in the right panel 



FIG. 21: A distribution of sky observation for an observa- 
tional strategy mimicking that of the EBEX balloon-borne 
experiment. The left panel shows the full observed sky patch, 
while the right one - its well-observed part. The latter con- 
tains pixels for which a number of the observations is not 
smaller than 5 x 10"*^ of the best observed pixel. The density 
of observations per pixel ranges from ~ 2 x 10^ in the center 
down to 1 or 10^ at the edge in the full and reduced sky cases 
respectively. The pixel size is roughly 7 arcmin, correspond- 
ing to Nside = 512. The size of the shown panels is roughly 
20° X 20° each. The color stretch adopted here is linear. 



parable to that of the Fisher variance at least down to 
i ~ 200. For larger angular scales however, the patch 
geometry starts affecting the performance of the estima- 
tion of the power spectrum as the variance increases for 
more elongated patches. A similar behavior can be seen 
also in the case of the optimized windows, as expected 
given that the patch geometry does not affect the level 
to which the E /B leakage is controlled. 



C. More realistic example 

The sky coverage of the realistic experiments is clearly 
more complex than any of the cases studied so far. The 
observed patches are usually irregular and a density of 
observations per sky area is pixel-dependent giving rise 
to significant noise inhomogeneity of sky maps recovered 
from data of such experiments. In this Section we there- 
fore assess the performance of the pure pseudo-power- 
spectrum approach in the case of an observation mimick- 
ing a long-duration CMB balloon-borne experiment, such 
as for instance EBEX [28]. The adopted here sky cover- 
age and its sampling, i.e., a number of observations per 
pixel, corresponds to one of the possible scanning strate- 



gies for that experiment and is shown in Fig. 21 As it is 
typical of the small-scale experiments the realistic distri- 
bution of the observations is very inhomogeneous. The 
number of samples per pixel ranges from 1 at the patch 
edge to roughly 2 x 10^ in the center. Though in princi- 
ple including all the pixels could be advantageous, if they 
are properly apodized, in practice it turns out not to be 
always the case. Moreover, we find that due the dynam- 
ical range for the noise levels the PCG-based, iterative 
calculation of the optimized windows fails to converge 
within a reasonable number of iterations, while a direct 



of Fig. 21 It contains all the pixels for which a number 
of observations is not smaller than 10^, reducing the dy- 
namical range of the noise levels from roughly 3000 to 
15. This has an unavoidable and unfortunate impact on 
a observed sky area, which in the studied case decreases 
from 1.1% down to 0.89% of the total sky. We set the 
noise level per sample in such a way that average noise 
level for the full patch corresponds to that used in the 
homogeneous noise cases studied earlier. 



In Fig. 22 we show the variances of the B-mode spectra 
computed using full (left panel) and only well-observed 
part of the patch (right panel). The gain with respect 
to the noise variance largely compensates for the loss in 
terms of sample variance due to removal of the very noisy 
pixels close to the edge of the full survey area. This is 
especially apparent for ^ < 150. When analyzing the 
full sky survey, long apodization lengths are preferred 
as underlined by the red crosses displaying the minimal 
variance achieved for each bin. This is because a long 
apodization length allows us better to mitigate the in- 
fluence of the noisiest pixels at the edge and thus to 
lower the noise variance. This also explains why small 
apodization lengths are suboptimal for the entire range 
of considered ^-modes. However, once the noisiest pix- 
els are removed (top-right panel), the situation changes 
drastically. Moderate apodizations are preferred at large 
angular scales (4 degrees for £ < 100). This value of the 
optimized apodization length is the same as for the ho- 
mogeneous case. Nevertheless, at small angular scales, 
the behavior of the variance with respect to the apodiza- 
tion length for the inhomogeneous noise differs from the 
homogeneous noise case. For the latter, small, ~ 2 de- 
grees, apodizations are preferred at small angular scales 
whereas long apodizations should be used at small scales, 
~ 8 degrees for £ > 800, for the former. Such a modifi- 
cation is due to the noise variance (we remind that with 
analytic windows, the E-to-B leakage is sufficently low- 
ered below the noise level). In this part of the spectrum, 
windowing should converge towards inverse-noise weight- 
ing for the estimation to be close to optimal. In the ho- 
mogeneous noise case, this requires to lower as much as 
possible the apodization length. In the inhomogeneous 
case, however, this means that windowing should be as 
close as possible to the actual noise distribution. For our 
realistic example, the noise distribution as shown on Fig. 
[21] exhibits a relatively high apodization length, which 
translates into a large apodization needed for the window 
to be optimal. 

In a more general case the question if any given ana- 
lytic window can produce the results close to the optimal 
is more complex and will depend on both the window 
and the noise distribution. In fact, at high ^, the near 
optimality can be reached only if the adopted window 
function approximates the inverse noise weighting. 
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FIG. 22: Upper panels: Variance of the estimated B-mode power spectrum computed using C^-windows with different apodiza- 
tion length (color curves) for the full observed (left panel) and the well-observed (right panel) sky area. The black-dashed 
curves show the Fisher estimate of the variance calculated assuming the average noise level (see text) and are only shown as 
a reference here. The black solid curve represents the input B-mode power spectrum. The red crosses indicate the lowest 
variance found for each bin, thus defining the optimized apodization length selected via the MC optimization procedure. For a 
given apodization length, removing the noisiest pixel significantly reduces the noise. Nevertheless, except for the first two bins 
the optimized variance is nearly the same for the two considered sky patches. Lower panels: Same as above but using a window 
function given as a product of the analytic C^-window and the inverse square-noise variance. This extra noise weighting largely 
alleviates the problem of the noisiest pixels resulting in similar variance estimates for the full and reduced observed sky patches. 
Moreover the optimized apodization length, marked with red crosses, is independent of the bin number. 



Our discussion and results derived above show that the 
analytic windows introduced in this paper can indeed be 
appropriate for noise patterns similar to EBEX. A more 
general approach, applicable to other noise inhomogene- 
ity patterns is also however possible. We first note that 
multiplying any of the analytic apodization studied here 
by some smooth function of the sky position will result in 
a proper apodization, which suppresses the E/B leakage. 
We emphasize that this observation is indeed strictly true 
only for nonpixelized skies and as usual a care has to be 
exercised whenever pixelization is introduced. In gen- 
eral, for pixelized skies the smooth function needs also to 
be sufficiently flat at the boundary not to lead to some 
uncontrolled level of leakage. In the case at hand we de- 
rive the new spin-0 window by multipling the analytic 
windows by a pixel-dependent weight corresponding to 
the inverse squared noise weighting, cx l/cr^(p), where, 
cr^ {p) is the estimated noise variance at pixel p. We then 



numerically compute the spin-1 and spin-2 components 
from it. This allows us to mimic inverse- noise weight- 
ing while fulfllling the proper boundary constraints. The 
effect of such inverse-noise weighting is displayed in the 
lower panels of Fig. 
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for the full survey (left panel) and 
its well-observed part (right panel) , where colors refer to 
the value of the apodization length of Wq. The results 
show that using this inverse-noise weighting signiflcanlty 
reduces the total variance, particularly for small apodiza- 
tion length. This explains the origin of the extra- variance 



at low-^ end seen in the upper panels of Fig. 22 relating 
it to the noise of the noisy pixels present at the map 
boundary. 

Interestingly whatever survey we use (full or reduced) , 
the lowest variance is always reached for a small apodiza- 
tion length {Sc = 2 degrees) apparently alleviating the 
need of going through the optimization procedure of the 
apodization length. We note however that this fact is 
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FIG. 23: Variance of the estimated B-mode cross-spectra using the C^- and C^- windows as well as their inverse-noise weighted 
extensions for the full observed (left panel) and well-observed sky (right panel) areas. The black-dashed curves show the Fisher 
estimate of the variance derived for the average noise level and the black solid curve depicts the input B-mode power spectrum. 
If the full observed map is to be analyzed, adding inverse-noise weighting suppresses the variance below the estimated power 
spectrum level at the low-^ regime. With the noisiest pixels removed, adding inverse-noise weighting does not change much 
except at large angular scale where it results in some excess variance, when compared with the simple analytic windowing. 



related to the specific noise distribution considered here 
and will be applicable in general only to the cases with a 
sufficient, noise-induced apodization present at the patch 
edges, which, as shown here, are however of a practical 
interest. At the high--^ part, this is because the near op- 
timality is reached thanks to the inverse-noise weighting 
and therefore the small apodization length of the analytic 
part of the window is preferred, as it affects the total win- 
dow shape to the least extent. At the low-^ part of the 
spectrum, without the extra noise-weighting an apodiza- 
tion length of 4 degrees was needed for the analytic 
weighting alone. The effective, noise-induced apodiza- 
tion is larger than that and the apodization length for 
the analytic part of the window can be therefore smaller 
as the latter is needed only to fulfill the boundary condi- 
tions and reduce the E/B leakage. 

Somewhat counter-intuitively the numerically derived 
estimates of the power spectrum uncertainties are seen 



in Fig. [22] to be superior to the forecasts based on the 
Fisher formalism. This is due to the fact that the latter 
are based on the average noise level computed assuming 
the same number of observed samples as in the actual 
case, but now homogeneously distributed over the entire 
observed sky area, and the numerical calculations bene- 
fit, thanks to either noise- weighting or appropriate ^-bin 
dependent apodization, from the very low noise core of 
the patch, which dominates the constraints at the high-£ 
end of the spectrum. 

A comparison of the best estimates of the spectral vari- 
ance derived with different analytic windows is shown in 
Fig. [23] It demonstrates the importance of the down- 
weighting of the very noisy pixels for the full patch anal- 
ysis (left panel), which is needed to achieve a variance 
level lower than the input B-mode at large angular scales. 
This effect is particularly prominent for the -window, 
which by construction is less smooth towards the patch 



boundary than the C^-window, and adding inverse noise 
weighting is here mandatory. If the noisiest pixels are 
removed from the analysis (right panel), using inverse- 
noise weighting performs the same as analytic windows 
alone except for £ < 200, where inverse- noise weighting 
is somewhat suboptimal as compared to analytic win- 
dowing alone. This effect is analogous to the usual loss 
of precision in the signal-dominated regime also seen in 
the standard pseudo power spectrum cases. However, 
as noted previously, the results derived with the explicit 
inverse-noise weighting have not required any MC-based 
optimization in the case studied here. Therefore though 
suboptimal for £ < 200 they can certainly be used to pro- 
vide a first and reliable estimate of the B-mode power 
spectrum, at least in the cases whenever the noise distri- 
bution complies with the assumptions as already spelled 
out earlier. 

In Fig. [24] we show the variances derived using the nu- 
merically optimized windows and compare them against 
those obtained with the help of the analytic ones dis- 
cussed already before. Only the reduced, well-observed 
sky area has been used here, due to the numerical con- 
vergence problems experienced by the PCG procedure 
as mentioned before. The apodization length of the an- 
alytic windows (which include no noise-weighting here) 
has been optimized in each bin in order to minimize the 
variance. Clearly, the windows optimized numerically in 
the pixel-domain display very similar levels of the uncer- 
tainty on the recovered power spectrum level as derived 
with the analytic apodizations. The harmonic space op- 
timized window shows however an excess of variance in 
the entire range of considered ^-modes. This is because it 
has been calculated using the average noise level, which 
in all bins has incorporated all the pixels, including those 
close to the patch boundary, therefore it does not benefit 
from the low- variance of the pixels in the very center of 
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FIG. 24: Variance of the estimated B-mode power spectrum 
for the reaUstic study case considered here. Four different win- 
dows have been used : C^-window (yellow curve), C^-window 
(red curve) , optimized window computed in the harmonic do- 
main (green curve) and the optimized window computed in 
the pixel-domain (blue curve). The apodization length of the 
analytic windows has been optimized in each bin. The black- 
dashed curve stands for the Fisher estimates of the variance 
and the black-solid one for the input theoretical B-mode. The 
Fisher estimate has been calculated assuming the average 
noise level. All the windows show comparable performance 
with an exception of the harmonic domain optimized one. 
This departure is due to the homogeneous approximation of 
the noise making the effective apodization length of the har- 
monic optimized window suboptimal. 



the patch. At the high-^ end the variance obtained with 
the harmonic domain window follows closely the Fisher 
matrix estimate, which has been derived using the same 
approximation for the adopted noise level. 



V. CONCLUSIONS 

We have implemented the pure pseudo-spectrum for- 
malism of [15, 16 extending it to incorporate the cross- 
spectrum approach and discussed practical issues related 
to its implementation. In particular, we have presented 
a quick, efficient way to compute pure multipoles and 
their respective mode mixing kernels utilizing harmonic 
domain derivatives, and demonstrated its overall consis- 
tency. We have also considered a number of apodization 
proposals. Those have included analytical, semianalyti- 
cal and numerically-derived functions. We have studied 
the performance of the implemented formalism from the 
point of view of a level of the residual E/B leakage as 
well as final variance of the estimated spectra, and inves- 
tigated their dependences on a number of factors such as 
an observed sky geometry, window choice, pixelization 
type and pixel size. We have considered the role of the 
sky priors in the apodization optimization procedures, as 
well as proposed an iterative, prior-free approach to the 
B-mode power spectrum estimation. 



Our results support the assertion of [151 [IS] 
show that the pure pseudo spectrum approach not only 
can successfully resolve the E/B leakage problem but, if 
properly applied, it also allows to bring down the total 
variance of the estimated power spectrum to levels com- 
parable with those derived from the optimal approaches. 
And it can do that at the significantly lower computa- 
tional costs. In fact, we usually find that the power spec- 
trum variance derived with the pure approach is within a 
factor of 2 away from the reference variance based on the 
Fisher matrix approach, which neglects the E/B leakage 
and assumes no apodization (but a binary mask). 

We have shown that the analytic apodizations, in par- 
ticular the (7^ -window proposed here, can deliver a com- 
petitive performance, particularly if optimized using the 
MC simulations, at the low calculational cost. The appli- 
cability of the analytic windows can be greatly extended 
by combining them with the inverse square-noise weight- 
ing. Such weighting may sometimes produce somewhat 
suboptimal results at the intermediate and low end of 
the ^-range, where it can/should be replaced by the un- 
weighted analytic or numerically optimized windows. We 
have however found that in some cases of practical inter- 
est the noise weighted analytic windows may not require 
any optimization procedure. They therefore provide a 
suitable first guess useful in particular in a preliminary 
quick-look analysis. 

In the cases with a nearly homogeneous or mildly inho- 
mogeneous noise distribution, we have derived a semiana- 
lytic, quick-to-compute harmonic domain window, which 
though requires some prior information about the sky sig- 
nal, provides yet another cheap and efficient alternative 
to be used in the pure spectrum calculation. However, in 
the highly inhomogeneous cases such a procedure may be 
suboptimal and further work is needed to bring its per- 
formance in-line with the other methods discussed here. 

The numerical window optimization as proposed in ^6] 
have been shown to provide an efficient framework for 
the optimized window computation applicable to a wide 
variety of the problems. It can however be numerically 
prohibitive and iterative solvers needed to overcome this 
problem may suffer due to lack of convergence. We have 
found that a particularly common case is whenever a high 
dynamic range of expected noise levels is present. In such 
circumstances a conservative tresholding procedure may 
have to be applied in practice to avoid such difficulties. 

We note that although these conclusions have been de- 
rived from the cross-spectrum analysis they should be 
equally applicable to autospectra cases. 

In this work we have not considered issues related to 
bin-bin correlations of the power spectrum estimates. 
This is an important topic for the further applications 
of the results of the method discussed here, in particular 
for the cosmological parameter estimation, and will be 
treated in the forthcoming work. 

This paper focuses only on the statistical uncertainty 
incurred due to either the noise and/or sampling vari- 
ance. Other potential sources of error, for instance, re- 
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suiting from a foreground component separation and/or 
systematic effect treatment are not included. Clearly, 
such effects will have to be incorporated in a final error 
budget of B-mode power spectrum estimated from any 
realistic CMB polarization data. The results presented 
here set a stage for an analysis of those issues in future 
work. 

In this work, we have not discussed the role and a po- 
tential impact of pixels missing from the main body of 
the observed map. Such "holes" in the maps could result 
from, for example, a point-sources removal or masking. 
We note that their presence does not lead to any fun- 
damental problems in the pure formalism, as has been 
shown in the case of the pixel-domain optimized windows 
in [16j. However, the practical issues exist. In particular 
in the cases with the analytic windowing, whenever the 
missing pixels are closer to each other and/or observed 
patch edge than the considered apodization length a spe- 
cial ca re must be taken to avoid windows discontinuities 
(Sect. IV A 4). Such effects may cause some additional 



E/B leakage, however they are unlikely to lead to any 
appreciable increase of the total variance of the recov- 
ered B-mode power spectrum. A more detailed study 
of the effect of holes in B-mode reconstruction will be 
presented in a future work. 

The results derived in this paper were obtained 
with the help of a publicly available software package 
S^HAT [30 . The set of MPI-parallel routines permitting 
quick computation of the pure multipoles as well as spin- 
weighted apodizations given an input spin-0 window, de- 
veloped in the course of this work, can be downloaded 
from ^37] . 
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APPENDIX A: MIXING KERNEL 



The mixing kernel is implicitly defined by Eq. (12). 



It can be explicitly derived in the harmonic space by 
expressing the pseudo-Q's, averaged over CMB realiza- 
tions, as functions of the CMB power spectrum. The 
derivation is very similar to the one described in [TT\. 
We first recall that the pure pseudo-a^^'s are given by 
the following expression: 



EO^im = E2/m + 2— + — Eo^£rn, (Al) 



Ne,i^ . 1 



BCi£m — ^2,^m + 2— ^^i^^^ -\- — — 
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where N£^s = \/(-^ + ^V-fi^ ~ ^V-- Employing now har- 
monic space decompositions for the polarization field 
Q -\- iU and the different spin-5 window functions, we 
can show that. 
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The above expressions can be then inserted into Eqs. (Al) and (A2), which we can use to define the pure pseudo-C^'s. 
Taking then the average over CMB realizations and making use of the orthogonality relations. 
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we arrive at, 



CP 



On defining tiie power spectrum of tiie spin-0 window function by 



Wf" = — - — V 

^ 2f ' + 1 ^ 

and introducing the following notation 

Jti^J J )= \ _2 + 5 2 -s I ^ \ 2-s -2 s 
we can express the diagonal and off-diagonal blocks (in the E/B subspaces) of the mixing matrix as: 
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We note that in the standard pseudo-spectrum approach, 
the mixing kernels are obtained by removing the and 
terms in the above expressions as the latter encode 
the contribution of the counterterms. 

If the window function is apodized up to its first deriva- 
tive, the off-diagonal block of the mixing matrix should 
be equal to zero. We have checked numerically that this 
is indeed the case, when the above formulae are applied. 

For the binned power spectrum estimation, the mixing 
matrix as implicitly defined in Eq. (12), is given by. 



In the above formulae we have assumed that the deriva- 
tive constraints between the three spin-s windows are 
exactly fulfilled. However, in practice they may be only 
approximately satisfied and the mixing kernel should be 
therefore expressed not only as a function of but also 
of Wi and W2. If such extensions are incorporated, the 
mixing kernel reads. 
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where 

w^srm"=\ W^s(S^)^sYln^,md'^- (A9) 

J47r 

These last expressions can be rewritten using the E and 
B decomposition of the spin- weighted window functions, 
defined as, 

^O^L -^0,£m, (AlO) 
I 



^s,tm (^5,^m + (-l)^^-s,^m) , (All) 



^i^m = ^ (^s,£m - (-l)^^-s,^m) , (A12) 

with 5 = 1 or 2. Using this decomposition in the expres- 
sion for the mixing kernel leads to. 
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We point out that because the spin-0 window is real, 
it does not have any B component. As a consequence, 
when the derivative relation is completely fulfilled, the 
spin-1 and spin-2 windows have also a vanishing B part 
and the second line in Eqs. (A13) and (A14) are equal 



to zero. This is the case when using analytic windows or 
optimized windows computed in the harmonic domain. 
However, because the derivative relation is relaxed dur- 
ing the pixel implementation of the optimized windows, 
the B part of the two spin-weighted windows is not zero 
anymore in this special case, though it remains much 
smaller than their E part. 

Equations (A13) and (A14) are directly applicable 
to the mixing kernel computations for both auto- and 
cross-spectra, but in the latter case only when the same 
apodization is adopted for all the maps. However two 
different apodizations are often required in the cross- 
spectrum case, for example,x because noise properties of 



the two maps differ, and then the above equations have 
to be appropriately generalized. This can be done by rep- 
resenting each of the moduli present on the right hand 
sides of Eqs. (A13) and (A14) as a product of a respec- 



tive expression, of which the modulus is taken, and its 
complex conjugate and subsequently using one of the two 
different apodizations for calculations of the expression 
and the other - its conjugate. For example, the square 
modulus 
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in Eq. ( |A13[ ), is replaced by 
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where 



and 



are the window func- 



tions apphed to the first and secon d map respe ctively. 
In this paper we use Eqs. (A13) and (A14), or their 



generalization as described above, for the calculations of 
all the mixing kernels. 



where C = (dd*) is the full covariance matrix of the data 
consisting of the two maps concatenated together, i.e.. 



-N, 



(11) 



APPENDIX B: OPTIMIZED WINDOW FOR 
CROSS-SPECTRA 

In this Appendix, we show how the optimized ansatz 
proposed in |16] in the context of the autospectra can be 
extended to incorporate the case of cross-spectra. De- 
pending on the general context, different choices can be 
adopted for such an extension. We will present three op- 
tions, focusing on the case of two maps, d*^^^ and d*^^^. In 
the following, i, j indices will run over pixels and A, B 
indices over maps. 

The first and simplest choice is to optimize with respect 
to the auto pseudo-spectrum of each of the maps. The 
optimized windows are then directly given by the ansatz 
of ^16j applied to each of the maps separately. 



No. 



(00)p(a)^(0) 



1, 



i=i 

Nobs 



(Bl) 
(B2) 



where C^-^ = (^dfd^^ is the covariance matrix of the 

A-th map and P^^^ the geometrical matrix projecting in 
the ^, m space (see [16 ). 

The second option optimizes with respect to the au- 
tospectrum of the two maps combined together in the 
optimal way. In this case the optimal, quadratic estima- 
tor and the pseudo-spectrum one can be written as. 



Nobs 1 A R 

E E (c-p'")c-)^. d 



ij = l A,B=0 

Nobs 1 

E E 



(B3) 
(B4) 



with S and N the signal and noise covariance matrices re- 
spectively. From the above definitions of the estimators, 
the formalism of [16^ is easily generalized to get. 



Nob 



1, 



(B5) 



-ij - ij J 

i,j=0 B=0 

This is easily generalized to an arbitrary amount of maps 
by summing over the entire set of maps. 

To compare the above optimized window to the one 
given by Eqs. (Bl) & (B2), we assume the noise proper- 



ties to be the same for the two maps. As a consequence, 
optimal windowing will be identical for each map. By 
plugging this into Eq. (B5), one can easily show that 



optimal windows coming from this last equation differs 
from optimal windows coming from Eqs. (Bl) & (B2) 
by the way the noise is accounted for: with Eq. (B5), 



the noise part is divided by a factor 2 as compared to 



Eqs. (Bl) & (B2). If Nmap maps are to be analyzed. 



then the noise contribution would be divided by a factor 



Nmap instead of 2. 



The third and final option we consider consists of 
optimizing the weighting with respect to each cross- 
spectrum. For two maps, the optimal, quadratic and 
pseudo-sepctrum estimators read. 



No, 



c 



£ df) (c-V")C-i)] dW, (B6) 

|j;df)w,(°^)p(fW,(°^)d«. (B7) 

(B8) 



The optimal, quadratic estimators can be recasted as a 
function of the inverse of the covariance matrix D = C~^, 



Nobs Nobs rj. 



(B9) 



i,j = l m,n=l 



The explicit expression of D are given by, 

q(00) _ q(01)q(11)-1q(01) 

-1 



2)(00) 



q(11) _ q(01)q(00)--^q(01) 

2)(oi) ^ _c;(oo)-1q(oi)2)(ii)^ 
2^(10) ^ _c(ii)"^c(oi)2^(oo)_ 



(BIO) 
(Bll) 

(B12) 
(B13) 



Following [16^ , the optimized window is obtained by min- 
imizing the distance between the optimal, quadratic es- 
timator and the cross pseudo-spectrum estimators with 
respect to W''^^\ This leads to the following linear sys- 
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tern, 



No 



EC 

2=1 

- y 

9 J 



(01)p(a)TT.(01) 
ik ^ik ^^i 



(01)t3C 



k [ 



D 



(00) 



D 



(01) 



(B14) 



This last approach is probably the more appropriate 
as the resulting windows are optimized for each cross- 
spectra. However, due to the rather complicated right 
hand side in Eq. (B14), it is also more numerically in- 



volved than the two first approaches. 



APPENDIX C: HARMONIC COMPUTATION OF 
OPTIMIZED WINDOWS 

In this Appendix, we describe the major steps to com- 
pute the optimized windows in the case of white and ho- 
mogeneous noise. As such noise is completely described 
by its power spectrum, we perform all the computations 
in the harmonic space. We start from Eqs. (18) & (19) 



and thus choose to use the x-fields rather than the spin- 
weighted approach in parlance of [16 . Consequently, we 
will optimize only the spin-0 window and compute the 
spin-1 and spin-2 windows as its numerical derivatives. 
We remind that the fi.e\-d is a scalar field defined as 

m 



(CI) 



where the ^-dependent prefactor compensates for the i- 
dependant normalization of the operator as defined 
in Eq. The harmonic decomposition of the fi.eld 
reads 



.B V- /(£ + 2)! 



(C2) 



where af^ denotes a type B multipole of the polarization 
vector P. 

Following the variational interpretation as proposed 
in [16 , the optimized windows are the one which min- 
imize the total aliased power in the pseudo-Q's. The B 
pseudo-spectrum averaged over CMB realization reads 
then, 



= Yl ^le-^oz-m-wl^.^n, (C3) 



where. 



_ + l(^-2)!(f + 2)! fie r 

~ ^ Att (£ + 2)!(f -2)! VO 

X {C^Bl^a'), (C4) 



and wo^irn stands for the harmonic representation of the 
spin-0 window function. This window has to be apodized 
up to its first derivative to make Eqs. (18) & (19) and 
Eqs. ([T6| and (17) (i.e., the x-fields and the spin-weighted 
approaches) equivalent. 



As a consequence, we have to minimize {Ci^ with re- 
spect to wo/rn under the three constraints. First of all, 
the spin-0 window has to be normalized within the ob- 
served region. 



y^UJo/ruMl^ = 1, 



(C5) 



£m 



where Aiim is the harmonic representation of the mask. 
(The normalization constant can be set to unity without 
loss of generality.) Second, the spin-0 window and its 
first derivative, the spin-1 window, have to vanish at the 
contour of the observed region, denoted V in the follow- 
ing. This translates into 2Nc constrains, where Nc is the 
number of pixels on the contour 



W(i) =0 and Wi(z) =0, for ah 2 G V. (C6) 



With such a large number of external constraints, the 
standard Lagrange multiplier technique used to derive 
the optimized spin-0 window function becomes hard to 
solve. The complexity of the problem can be reduced by 
assuming that the integral of the spin-0 and the spin-1 
windows on the contour have to vanish, leading to only 
two additional constraints 



Ewo/„^^L=0, (C7) 

VI{lTl)wo,£mVL = 0- (C8) 



£m 



Here Vim is the harmonic representation of the contour. 
The second equation above co rresponds to the spin-1 con- 
straint as wi^irn = V^i^ + l)'^o,£m- The optimized spin- 
window can be then derived by applying the standard 
Lagrange multiplier techniques, resulting in. 



mi" 



(C9) 



The two constants, A and /i, are derived by plugging the 
above solution into the constraints. Defining the follow- 
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ing total power, 
Na,o = 

Na,i 

Na,2 

Nx,o 
Nx,i 

Nx,2 



7^ \Me"^"\ 



e"m' 



X] 7^ \'Pe"m"f , 

'^e on 



^ P on 



m" I 9 



i"m" 

E 

l"m" 

E 





-1) 






-1) 



i"m" 

the two constants read, 



^ oil 



l2 

'm" I •) 



(CIO) 
(Cll) 
(C12) 
(C13) 
(C14) 
(C15) 



(C16) 



A,U^A,2 



Na,iNa,2 - N%^^ 

^ [{NxfiNAfl - Nx,iNx,2) 

,2, 



yW+^{Nx,oNx,2 - Nx,iNa,i) 



(C17) 



In this derivation the boundary constraints have been 
partiahy relaxed, making the window only partially op- 
timize. The numerical experiments, however, show that 



the resulting power spectra uncertainties are at the same 
level as the ones obtained using the exact, pixel-domain 
computation of the optimized window (see Fig. 18 ). This 
demonstrates that the proposed boundary constraint 
simplification does not compromise the performances of 
this fast computed, optimized window. 

We point out that the above approach is general and 
does not apply only to the 5-mode fields. For instance, 
we can use it to calculate optimized windows for a tem- 
perature (scalar) map. It suffices then to replace 



(^-2)!(f + 2)! 
(^ + 2)!(f -2)! 



in Eq. ( C4 ) by the total (signal plus noise) power spectra, 
i.e., CiiB'f, -\- a'^, of the map. Moreover, we can straight- 
forwardly relax some of the boundary requirements. For 
example, if no boundary conditions are to be imposed, 
we just need to replace Eqs. (C16)-(C17) with. 



A 



0, 



A,0 



(C18) 
(C19) 



while in a case when only the Dirichlet boundary condi- 
tion is to be imposed, the two constants are given by: 



A = 



Na,i 
Na,oNa,i - N^^o 

Nxfi 
N^o-Na,oNa,i' 



(C20) 
(C21) 
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